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PERTURBATION SOLUTIONS OF THE 
ELLIPSOIDAL WAVE EQUATION 


By F. M. ARSCOTT (Makerere, Uganda) 
[Received 1 November 1955] 
1. Introduction 


WHEN we separate the wave equation in ellipsoidal coordinates, we are 
led to the following ordinary linear differential equation: 


d*w 


a (a+bk*sn?z-+-qk'sn4z)w, (1) 
dz? 


where the Jacobian elliptic functions are formed with modulus k; q is a 
determinate constant, depending on the wave number, and a and 6 are 
arbitrary constants which we have to dispose of in such a way that the 
equation (1) shall possess a doubly-periodic uniform solution with real 
period 2K or 4K, and imaginary period 2:K’ or 41K’. It is to such 
solutions only that I shall give the name of ‘ellipsoidal wave functions’; 
solutions with Other periods exist but do not normally satisfy the 


necessary conditions for being physically significant. 

The equation (1) has been known in this context certainly since 
Méglich’s paper of 1927 (1) but has proved considerably less tractable 
than other equations with periodic or doubly-periodic coefficients, and 
the only other substantial contribution to the theory of such functions 
has been made by Malurkar (2). Apart from a short section in Strutt’s 
‘Lamésche, Mathieusche und verwandte Funktionen’ (3) and one paper 
by Campbell (4), no further work appears to have been done on the 
equation. It is, however, of considerable interest from the purely mathe- 
matical point of view, apart from its physical significance, since it 
includes, as special or limiting cases, such well-known equations as those 
of Lamé and Mathieu, the associated Mathieu equation (or spheroidal 
wave equation), and Hill’s equation with three terms. 

The perturbation method of solution, historically the earliest used 
for the solution of differential equations with periodic coefficients, goes 
back to Mathieu’s work in 1868 on the equation which bears his name. 
It consists in assuming a solution of the equation, possessing the neces- 
sary periodic properties, in the form of a series of powers of one of the 


Quart. J. Math. Oxford (2), 7 (1956), 161-74. 
3695.2.7 M 





162 F. M. ARSCOTT 


parameters of the equation (in this case q), the coefficients of this series 
being functions of the independent variable; at the same time the 
other parameters of the equation (in this case a and b) are also assumed 
to be expressed as power series in the first parameter g. Having made 
these assumptions, the coefficients in the various series are found suc- 
cessively by rejecting all terms not possessing the necessary periodic 
properties. 

The essential disadvantage of the perturbation method is that it is 
very difficult to determine the convergence of the resulting series. Even 
for the comparatively simple perturbation solutions of Mathieu’s equa- 
tion the problem has not been fully solved, and it is not yet possible to 
state any lower limit to the radius of convergence of the series which will 
be obtained in this paper. For this reason, the use of perturbation 
solutions for the other differential equations with periodic coefficients 
has now been largely abandoned, particularly with the development of 
the technique for dealing with series whose coefficients are given by a 
three-term recursion system. In the case of the equation (1), however, 
no such series has yet been discovered, and, although the author in some 
unpublished work has developed a somewhat similar series, it makes 
extensive use of Lamé polynomials which so far are very inadequately 
tabulated. 

In spite, therefore, of the disadvantages of the perturbation solutions, 
they are, in fact, virtually the only explicit solutions of equation (1) so 
far known. In this paper, I give not only the results so far obtained, but 
sufficient of the intermediate working to make further development of 
the solutions as easy as possible. 

In the heavy calculations which this method involves, the principal 
feature is the indefinite integral | sn*"z dz for positive integral values 
of n. From the formula 


d 
P {sn”-!zenzdnz} = (n—1)sn"-*z—n(1+-k*)sn"z+ (n+ 1)k*sn"+2z, 
dz 


and writing | sn*"z dz = I,,,, we obtain the recurrence formula 
(n—1)I,_.—n(1+k*)Z, + (n+ 1)k*2,,. = sn”™-!zenzdnz, 
by which we can express J,, 4, etc., in terms of J, and J,. In fact, 


z—E 
h = 3", 
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but, since these terms are not doubly-periodic, their precise values are 
not needed. 
Using (2) we obtain the following expressions, where for brevity we 
write s, c, d for snz, enz, dnz. 
3k], = sed+-2(1+-k*)I,—J, (3a) 
4(1 1 k?) 8+7k?+-8k4 4(1+-k?) 
‘. sed 4 ———— ———— J, 3b 
3h? v3 32 =~ (3) 


bo _ .  6(1+k) 24-4 2324 24k 
Tk? I, s°cd + Ae ee 
8(1-+k2)(6—k2+6k4), 244 23k24- 24h4 
a re 15k# si 
1+k* 48+ 47k?+ 48k4 
Si ) sed +‘ — a sed + 
_| p-2 >; __ -2_1. 16]-4 
4, MLE R16 — K+ 16K4) 
35k ' 
+99k4+ 1046+ 12848 I 
35k6 ’ 
4(1-+-k2)(16—k2+ 16k) 
eae 35k6 
+ 79k2-+- 80k4 
80-+- 79k?-+- 80k s8cd-+ 


5k? I, skcd I, — 


sed +- 





‘i 8 
s'cd + 





128+ 104k? 








: 

: 10(1+-k?) . 
s®ced +- cle s’ed + 63k 
— 16(1-+k2)(30—k2+ 30k) 
7” 315k6 


640-4 592k*+ 579k4+ 592k*+ 640K9 





seed +- 





d +- 
315k 





(1+ k?)(1280— 256k? +- 1206k4— 256k* + 1280k8) 
—- I,— 
315k8 





a oe 


640+ 592k2+- 579k4+ 592k*+ 640K8 
oF 315k8 ty. its 





We notice (it is useful later) that the coefficients of scd and J, are in every 
case the same but with opposite signs. 


2. Notation 
When g = 0, the equation (1) reduces to 
d*w 


> (a+ bk*sn2z)w, 
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which is Lamé’s equation; it is known that, for this equation to possess 
doubly-periodic uniform solutions with periods 2K or 4K, 21K’ or 41K’, 
we must have b = n(n+-1), where n is an integer which may be taken 
positive, and a must have a characteristic value, of which there are 
2n-+-1 for a given value of n. Then the required solutions are the well- 
known Lamé polynomials; the solutions of (1) which we require must 
therefore be such as to reduce to Lamé polynomials when ¢ = 0. 

It is possible to prove that (i) to every Lamé polynomial there corre- 
sponds one and only one solution of equation (1) which is doubly-periodic 
(with periods 2K or 4K, 2iK’ or 41K’) and uniform, and which is such 
as to reduce to that Lamé polynomial when q = 0, (ii) there are no other 
solutions of (1) which are doubly-periodie (with these periods) and 
uniform. In other words, the correspondence between Lamé polynomials 
and ellipsoidal wave functions (as we have defined them) is one-to-one. 

(The proofs of these statements have been obtained by the author in 
some so far unpublished work, but are lengthy since they depend on the 
construction of ellipsoidal wave functions by a quite different method 
from that used in this paper.) 

This makes it possible to base our classification of, and notation for, 
ellipsoidal wave functions on the classification and notation used for 
Lamé polynomials; it is convenient, however, to make first a slight 
extension to the standard notation for the latter. In proposing this, the 
notation used by E. L. Ince (5) has not been overlooked; the latter, 
however, was introduced primarily to deal with solutions of Lamé’s 
equation when » is not integral, and, while applicable to the Lamé poly- 
nomials, is somewhat cumbersome when so applied. The predominance 
of the polynomials, both in practical applications and in the theory of 
ellipsoidal wave functions, is so great as to justify, in the author’s view, 
a separate notation. In any case the correspondence between Ince’s 
notation and that used here is quite simple and is shown in Table II 
below. 

By Lamé polynomials are meant those solutions of 


Iw 
a ii {—h+-n(n+ 1)k*sn2zhw (4) 


dz? 


which are of the form sn’z cn’z dn‘'z F(sn*z) where r, s, t may be 0 or 1 
and F(snz) is a polynomial in sn*z. There are thus eight different types, 
distinguished by the values of r, s, t and consequently by properties of 
periodicity and parity. I shall retain the general symbol £”(z) to denote 
any Lamé polynomial of degree n, but distinguish the eight types by 
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prefixing, when necessary, the letters u, s, c, d, sc, sd, cd, sed according 
as the function F(sn®z) is multiplied by unity, snz, enz, dnz, snzenz, 
snzdnz, enzdnz,snzenzdn z, respectively. For example, an expression 
of the above form with r = 0, s = t = 1, would be written as cd H(z). 
! specify the upper index m as the number of zeros in the range 
0<2z< K. 
It is often more convenient, however, to write the lower suffix, not as 
n but as 2n, 2n+-1, 2n+2, 2n+3 according as the polynomial is of the 
first, second, third, or fourth species. Then the eight types, with their 
elementary properties, are denoted thus: 


TABLE I 
Type Parity Periods Type Parity Periods 
ubm Even 2K, 2iK’ scE},+2 Odd 2K, 41K’ 
she. Odd 4K, 2K’ sd EM . > Odd 4K, 41K’ 
cEm ., Even 4K, 41K’ cd EM. Even 4K, 21K’ 
dE.) Even 2K, 41K’ sed EX}. . Odd 2K, 21K’ 


The advantages of this numeration are that 
(i) the polynomial F(sn?z) is in every case of degree n in sn*z; 
(ii) theranges of then, mareinevery casen = 0,1, 2,...,.m = 0,1, 2,...,”5 
(iii) each of the above functions possesses m zeros in 0 < z < K and 
n—mzerosin K <z< K+iK’. The former statement is the defini- 
tion of m, the latter follows from Erdélyi’s results in (6). 


The correspondence between this notation and that of Ince is easily 
seen to be as follows. The upper index in Ince’s notation denotes the 
number of zeros in the range —K < z < K; constant factors are ignored. 


TABLE II 


Eom =uER; Ea, = dE&,,; Esi@+? = ecE™,,.; Esi:? = eed EY, 
Ec) = cdE®,9; Eci™:} = cE®,,; Esi®*+} = dE ..; Es3™?} = sE%,.,. 


With this notation for Lamé polynomials, we simply classify ellipsoidal 
wave functions according to the Lamé polynomials to which they reduce 
when gq = 0. I use the general symbol el(z) to denote any ellipsoidal 
wave function and el”'(z) to denote that ellipsoidal wave function which 
reduces to H™(z) when g = 0. I then further prefix the letters w, s, c, d, 
etc., to the symbol el?"(z), when needed, in exactly the same way, so that, 
for example, sd el?(z) is that function which reduces to sd £3(z) for g = 0. 
We thus have eight types of ellipsoidal wave function, whose properties 
of periodicity and parity are the same as the Lamé polynomials of the 
corresponding types. 
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The values of a and 6 corresponding to a given ellipsoidal wave function 
can be denoted (though it will often not be necessary) by using the same 
n and m suffixes and putting the letters denoting the type in brackets, 
thus for the function sd el?(z) we denote them by a®®2 and b&®%, 

The functions thus defined are determinate except for a constant 
factor. For the calculation of perturbation solutions (though not for the 
development of the general theory) the most convenient method of 
normalization is as follows. If #”(z) and el?"(z) are written in the forms 


N 
sn’z en*zdn'z ¥ A, sn®?z 

—" Pp 

p= 


x 
and sn’zen*zdn'z } B,sn*?z 

p=0 
respectively, then we normalize £7"(z) by the stipulation that Ay = 1, 
and el?"(z) by the stipulation that B, = 1 for all values of q. 

The calculation of the first ellipsoidal wave-function of each type, 
that is to say those which for g = 0 reduce to 1, snz, enz,..., snzenzdnz, 
respectively, is simpler than the more general case because (as will be 
seen) it involves only integration while the more general case involves 
finding particular integrals of certain linear differential equations. We 
therefore deal first with this special case. 


3. Perturbation solutions for the first ellipsoidal wave function 
of each type 
Consider first the construction of wel®(z). For q = 0, this reduces to 
E%(z) = 1, and a and 6b both reduce to 0. We therefore set 


w=1+) ¢f(z), a=) a¢, b= > d¢, 


and substitute in (1). Comparing the coefficients of q, q?, q°,..., we obtain 


the sequence of equations 
f(z) = ¢, +d, k*sn*z+k* sn‘z, 
= (c,+d, k* sn*z+-k' sn4z) f,(z)+-c,+d, k* sn*z 
= f(z) fi(z)+¢.+d, k* sn*z, 
f3(z) = (e, +d, k? sn?z+ k4 sn4z) f(z) + (co+d, k? sn®z) f,(z)-+c,-+d, k? sn*z 
= folz)fil2) +fil2{fo(2)fil2)f 1(2)} +¢3 +s k* sz, ete. 
We solve these successively using the formulae (3), and choosing the c’s 
and d’s in such a way as to exclude from f(z) all terms which are not 
doubly periodic. From the first equation we have, using (3a), 


filz) = (4. —4*) + fd, + 3(1+%*)} 4+ 44 snzenzdnz-+constant; 
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but, if f,(z) is doubly-periodic, f}(z) must be doubly-periodic also and 
the coefficients of J, and J, (which are not doubly-periodic) must vanish. 
The constant of integration must also be zero or f,(z) would not be 
doubly-periodic but would contain a term in z. Hence we have 

c, = 1, d, = —3(1+k), = f(z) = 4k’ snzenzdnz. 
Integration gives Si(z) = $k sn*z+ constant, 
but the constant of integration must again be taken as zero or the 
normalization condition will be broken (in fact, none of the f(z) will 
contain a constant term). 

The values of f,(z), c,, d, thus obtained are then substituted in the 
next equation and f,(z), C2, d, similarly calculated. There are no note- 
worthy features of the calculation, and the results are given in full in 
Table III. 

To apply the same process to the construction of the first function of 
each of the other types, it is convenient to make first a preliminary 
transformation of the equation to remove the extraneous factors sn z, 
enz, dnz. These transformations are as follows; w represents a solution 
of the standard equation 


d*w : 
| a2 25_) 4 4. 
> (a+ bk? sn2z+-qk* sntz)e. 


If w = ysnz, then 


d (sute 4) = [(1+k2+a)sn2z-+ (b—2)k? sn*z+-gk4 sn®z]y. 


dz\ ~ dz 


If w = yenz, then 


d ie = 
dz\ ~ dz 


- [l+a+f{(b—2)k?—1—a}sn*z+- (2—b+-gk*)k? sn*z —qk* sn®z ly. 
i j q q y 


i 


If w = ydnz, then 


d dn2: “4 
Al "ke 


/ 


= [a+k*+ (6—2—a—k*)k* sn2z+ (q—b+-2)k* sn*z—qk* sn®z]y. 


If w = ysnzenz, then 
; (sn*s cen2z dy 
dz ae a dz 
= [(4+k?+a)sn2z+{(b—6)k?—4—k*?—a}sn*z+ 


+ (qk?—b+ 6)k* sn®z—qk* sn8z]y. 
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If w = ysnzdnz, then 
- [ant dn?z A 
= [(1+ 4k*+a)sn*z+ (b—7—4k*—a)k? sn*z+ 
+(q—6+6)k* sn*z—qk*sn®z|ly. (Se) 
If w = yenzdnz, then 
¢ (ent: dn*z = 
dz dz 
= [(1+/?+a)+{(b—6)—(1+-)(1+k?+<a)}sn*z+ 
+{qk*—(b—6)(1+k*)+1+k?+a}k? sn*z+ 
+{b—6—q(1+k*)}k* sn®z+-qk* snéz]y. (5f) 


d ae ee 
—|{sn*z en*z dn?z — 
dz dz 
= [(4+-4k?+-a)sn*z+{(b—12)k®?—(4+-4k?+a)(1+-k*)}sn*z+ 
+{gk?—(b—12)(1+k*)+4+4k?+a}k? sn®z+ 
+{b—12—q(1+k*)}k4* sn§z+-gk* sn™z]y. (5g) 
The calculation follows in every case the same general lines as for the 
function u el§(z), but there are certain devices which are not obvious and 
which shorten the working considerably; I shall therefore illustrate them 
by working through the expansion of d el?(z) in some detail. 
Since dE}(z) = dnz, and since, for dnz to be a solution of (1) with 
q = 0, we must have a = —k* and b = 2, we assume the solution of 


(5c) in the form 

y=1+)> ¢f, a= —kP+ 3 a¢, b= 2+ 3 4,¢, 
substitute in the equation and compare coefficients of powers of g. From 
the coefficients of q we have 


4 (dn?z f) = c,+(d,—c,)k? sn?z+(1—d,)k* sn‘*z—K sn®z, 
dz 
i.e. dn*zf, = ¢,+(d,—¢,)k7,4+ (1—d,)MI,—k*s,. (6) 

Instead of finding c, and d, first, we assume f, = Asn*z, where A is 
a constant to be determined, so that 


dn*zf = 2Ascd—2k*Astcd. 
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Now we expand the right-hand side of (6), using formulae (3 a), (3 b), and 
compare coefficients of s*cd, sed, I,, and J,, getting the four equations 
—2Ak? = —}k4, 
2A = }(1—d,)k*—£h(1+-F), 

0 = (d,—c,)k®+$k*(1+k*)(1—d,)—},°(8+ 7k + 8k), 

0 = c,—}(1—d,)k*#+4F(1+FP). 
The first of these gives A = ;jk*; the second then givesd, = —32(1+2k*); 
c, could be obtained from either of the last two equations, but much 
work is saved by combining the second and fourth equations, giving 
c, = 2A = 1k*. The third equation may then be used for checking. The 
possibility of eliminating d, between the’second and fourth equations 
results from the fact noted that the coefficients of J, and scd in the 
formulae (3) are numerically equal. 

(This procedure is possible only when expanding functions which do 
not contain snz as a factor. It appears to be true in the other cases 
that c, = 6A, but this is not immediately apparent.) 

Having thus found the values of c, and d,, and that f, = 4k? sn*z, we 
use the next equation, which is 

pre ~—re ae $a 
q, (an*zf 2) = C,+(d,—c,)k* sn*z—d, k*sn ethics (dn zf x). 
To solve this, we assume 
fo = Bsn*z+Csn*z, 
and on substituting this and the known expression for f, we compare 
coefficients of s'cd and s°ced, giving 





6 este @Lae 
—4i2B — o. 2(2B—k2C) = ka( 6( 4 k yf 4k ) 
(0 \ 350 250 
These easily yield 
4 2/2 1 ].2 
‘a 2” 8 ee 


280’ a 
From the two equations arising from the coefficients of scd and J, we 
at once find je 2 
.2(3 1 [2 
Co = 2C a A tt att I 
4 875 
Then either of these two equations may be used to find 


__ 9(2__ 2]).2_1 D]}-4 
iim 2(3—3k 2kA) 





875 
The remaining equation, that obtained from the coefficients of J,, may 
be used for checking. 
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Proceeding on the above lines, the following eight solutions have been 
developed as far as the term which is O(q?), and in the case of wel? as 
far as O(q°). There seems to be no difficulty apart from the complicated 
working in extending them further. The work becomes more compli- 
cated with the greater number of extraneous factors involved; for 
example, in calculating the term which is O(q?) in wel§ we need only the 
formulae (3a) and (3b): that is to say, the formulae for J, and J,; to 
calculate the corresponding term in scd el we need the formulae up to 
and including J)». 


TaBLe IIT 
[The terms omitted are O(g*) in every case except (1) when they are O(q*)] 


gk? sn@z k’‘sntz  k*(1-+4 | 


,]0/ + ! en 
(1) welg(z) = 1+—— +9" Gag 4 270 





_ gf Afsn®z | k4(1+k*)sn*z | 
T?'|20°39.5.7 22.3°.5.7 3.5.7 


a > 


wo — oH, RMI +E) 241 —4k* +284) 
bes Ist? a antl > 


PBF 





poe — —g 2+) 924K) 21 +k*(2— 5k? + 2k4) 
i ia eet ieee 3°.5.7 . 

= sn* sn‘z Px l | 
280 1750 ‘ 


4 
(2) self(z) = snz+q— = 4+-@@snz IF 


3k? 3k 1 +k) 


a®? = —1—k?+q—+ — 
3 1 q 875 


Lp D 2) OL 
bie — gg H(i +h)_ kt) 





k* sn@z ~ fa k*(1+ te 
10 +o _—”l—lC( iar 
ie 2, k%™(1+3k) 


875 


(3) cel}(z) = enz+qenz 


(c)0 __ 
as) = 





pee — 2 


k? sn2z kts _ M(8+k*)sn% 
ee ne ee BE pencils 
(4) del{(z) = dnz+qdnz 10 e144 dn - 1750 |. 


.2 2(3 4 2 


’ 2(1+ 2k?) 2(3— 3k?+ 2k*) 
(d)o _. 9__ __ 
ie dl it eee q 875 ; 


> 
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TABLE III (cont.) 


k? sn*z k'sn4z . k?(1+3k?)sn2z 
snzenz+ qsnzcnz—— --+-g*snzenz . 





50a 2.3.78 
_. Sk k(14-3k%) 
4 en ne eh el cen, 
B+ te — saa? 
2(2k24+3) A 1—k* +. 2k4) 
— a". 
3.73 


, ke? sn*z | e , i= sn*z | 
snzdnz+qsnzdnz——— +q’?snzdnz “> 
: ie be ey a i Os 


6—q 





k(3+k?) 
= te — 


a7? 
2(2—k?+ k*) 
k* sn®z 


k 
enzdnz-+4 qenzdnz—- genzdna| 


F e1 1 k?) 
a 


(2—3k?+ 2k*) 
3.73 





4snfz k2(1+ | 
504 + 2.3.7° ; 


-¢ 








> 


k? sn®z 


sed el§(z) = snzenzdnz+qsnzenzdnz “—_ 


k*sn*z k*(1+ | 
23. 3?.11 2.3°.11 P 





+ gsnzenzdn2| 
ae 2 k*(1+k*) 
7° vu 
2 2(1—k*+ k*) 
3*.11 





4. Perturbation solutions for other functions 

As has been mentioned above, the perturbation solutions for ellipsoidal 
wave functions other than those in Table III are more complex. The 
complexity arises, however, less from the solution itself than from the 
fact that the Lamé polynomials themselves are complicated expressions; 
this makes the working very difficult in the general case, where the value 
of k? is arbitrary. If, however (as is normally the case in practical 
applications) the numerical value of k? is known a priori, much of the 
difficulty disappears. 

I shall illustrate the calculation by considering the function w el}(z). 
The Lamé polynomial to which this reduces for q = 0, namely wH}(z) is 
1+#+,/(1-+K) (7) 

3k? 





uES(z) = sn®*z— 
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which satisfies the equation 
d*w 
dz2 


For brevity, write 


[—2(1+k?)+2,/(1—k2+i4)+ 6k? sn®z Jw. 


, is 1+k®+A 
A= ,(1—F+2), —G= me 


so that wH%(z) = sn?z—Q, satisfying the equation 


dy = —2(1+k*)+22, 


tn = (a, + 6k? sn*z)w. 
It is easily seen that 
32Q2—2(1+k2)Q+1=0, aQ= 
We assume the formal solution 
welf(z) = uEY(z)+ ¥ ¢f(z), 
a=a+>o,¢, b=6+)34,¢, 
and substitute in (1). Reducing and comparing coefficients of powers of 
q yields the equations 
fi—(ao+ 6k? sn2z)f, = (ce, +d, k* sn?z+e sntz)u EY(z), (9) 
f—(@o+ 6k? sn2z) f, = (c,+d, k* sn?z+ k* sn4z) f, +(c.+d, k* sn?z)uEY(z), 
etc. (10) 
When the right-hand side is expanded, (9) becomes 
fi—(ao+6k? sn2z) f, 
= —¢c, Y+(e,—d, k*Q)sn2z+ k*(d, —k*Q)sn4z+k*sn®z. (11) 
Now the normalization condition implies that, in this case, the co- 
efficient of sn?z in wel§(z) must be unity for all values of g, and hence 
fi. fo, ete. will not contain any term in sn*z. We therefore assume a 
particular integral of (11) in the form f, = Asn‘z+ B. Substituting and 
comparing coefficients of powers of snz, we have 
14k24 = kA, fay +16(1+k*)}4 = k*(k®?Q—d,), 
12A—6k°B = c,—d, k*Q, a,B = ¢,Q. 
Solving these in succession and substituting for Q and ay, we have 
A= iF, d, = —§(1+F*)+3, y= ee. 
(1-+-k2)(— 10+ 16k2— 10k4) +-A(— 10+ k2— 10k4) 
" 2.39. 7h 





B 
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Thus we have a particular integral of (9). The complementary function 
is aB3(z)+BF9(z), where F(z) is the standard second solution of Lamé’s 
equation corresponding to E£%(z). Since f, is to be doubly-periodic, and 
also to contain no term in sn®z, we reject this complementary function 
completely and have f, = Asn‘z+B, where A and B have the values 
given above. 

The expressions already obtained for f,, c,, d, may now be substituted 
in (10) and the same process continued. The work is obviously very com- 
plex, and it may be doubted whether the practical value of general 
formulae are worth the very heavy algebraic manipulation involved. 
If, however, the numerical value of k? is known, the work becomes 
comparatively simple since A, Q, a, can at once be given numerical values 
also. 

It is worth noting, however, that wel}(z), a™}, b™3, can be obtained 
from uel$(z), a\8, b§, merely by changing the sign of A. 


Checking 

The solutions obtained above may, of course, be fully checked by 
substituting the series for el(z), a, b in the equation (1). This, however, 
is a lengthy process, and a much simpler check is provided, at least for 
the expressions for a and b, by means of certain transformation formulae 
obtained by the author on lines similar to Erdélyi’s work on the corre- 
sponding formulae for Lamé functions [(6) § 5]. 

If, in the expressions for a, b, we write —q for q and k’? for k* (where 
k2-4-k’2 = 1), I shall denote the results by a, b, respectively. Then the 
transformation formulae are 


F(uym (ujn—m __ « =(u)m (umn Lg_(ujn—m __ 9 
Ben — ben, is | Gn +§ T43n => (12a) 


},(s.d)m (d,s ym—m __ 9 q@is-dm 1 f(s.dm (d,s )n—m __ 9 
OST Ost , Aon+1 +b¢,0"+ ald. +1 me Gs, (12b) 


Jlcym (c)n—m ‘ =(c)m (c)m i gle)n—m __ 9 
OORT - Aan+1 +05 iT? = ©, (12¢) 


bfec-edym __ 5 fed.seyn = 2g, aseee ,cd)m _| +-bi(see em | +a: . m — = ¢, (1 2d 


> 


b{edm _ p{ed)n— Ce aedm + | b(ed)m 1 g(ed)n— _ (12 e) 


},(sed)m__ ,(sed)n—m __ 9 ete bieod)m foot n— —m __ 9 
Ons —Oan +3 = <¥ ans 3 +3 +4sn+3 = (12f 


In (12b) and (124d), either the first, or the second, character must be 
taken throughout. As an example of the use of these formulae, let us 
take (12b) and put n = m = 0. Taking the first character gives 


Be _YO2 — 2g, Ges Hons gia — g, (13) 
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Putting —q for g and 1—k? for k? in the formulae (2) of Table III, we 


have , , 
3—3;° 6—9k?+-3h4 


a(so __ )_1 }2 
a} = —2+k —— 
875 


- __ 4p2 2(2__ 9,2.) op 
¢q8=tt ) 4 23 3k2+4-2k4) | 


bee a 2 
5 875 





and it is now a simple matter to verify that both the equations (13) are 
satisfied. 
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THE BEHAVIOUR OF INTEGRAL FUNCTIONS 
DETERMINED FROM THEIR TAYLOR SERIES 


By J. CLUNIE (Keele, North Staffs.) 


[Received 30 November 1955; in revised form 25 June 1956] 


1. In this paper I shall prove the following theorem. 


THEOREM 1. Let f(z) = Sa,2" be an integral function and « be a 
n+I1\*| a, | . ; i ik : 
—— fa | is ultimately a steadily increasing 
An+1 


positive number. If ( 
n 


function of n, then 
M(r,f) < {1+0(I)}a-2 AP (1 +aev(r,f alr, f )- (1) 
THEOREM 2, Jf f(z) satisfies the conditions of Theorem 1, then 
: log M(r,f) ' 
lim sup ————"~ < 1+ a". (2) 
pee? loguir.f) 


The above theorems are easily proved and their value depends on 
their being best;possible. This is shown by 


THEOREM 3. Given « > 0, there exists an integral function 


g,(z) = > On 2 


b 
M(r,g,) > {1—o(1)}a-®47T(1+-ae(r, 9)u(7, Ja) 
for an infinite sequence of r, and 


log M(r,f) ciao (4) 


n 


such that te ) >, | is a steadily increasing function of n and 


n+1 | 


lims ‘al 
oe. logulr.f) 


From the construction of g,(z) of Theorem 3 it will be seen how an 


integral function A(z) = > ¢, 2" can be found such that fn steadily 


n+1 
increases with n and 


lim up = 00 
mate {v(r, k)}Fu(r, h) ' 
log, M(r,h) _ 


lim sup —2=——-"—S = 00 
root pr, h) 


> 


Quart. J. Math. Oxford (2), 7 (1956), 175-82. 
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where k is any positive integer, however large. Hence there can be no 
general results corresponding to Theorems 1 and 2 for the case « = 0. 

At the end of the paper I shall use Theorem 2 to extend a result due 
to Singh and Shah (1). 


2. The proofs of Theorems | and 2 depend on an attenuated form of 
the Wiman-Valiron method [(3) 93]. If f(z) satisfies the conditions of 
Theorem 1, then 

— a 
$,(z) = S —B2" 
Pal?) aan 1 
1 
is an integral function for which all integers beyond some definite integer 
are central indices. Let v = v(r,¢,) and then 


Ay r? <— (a, ad (n = 1 2 ) 
~ 9 Syeee)s 
n= = y* 


ae a 
so that “a <— (n=1,2 
a, r’ yp* 


¥ 
: 

~— 
: 


For any v, R (< 1) can be chosen so that, for all n, n*R" < v*R’, and 


hence a,,\(rR)" _ n*R" _ : 1s) 
a,\(rR)Y ~~ v=RY ~~ ” ’ 
which gives v(r,d,) = v(R, F,) = vi(rR,f), (6) 
where F,(z) = s n*z", Further, 
T 
(rR, f) = \a,\(rR)’ 
= sl pw ya Re 
Vv 
= w(r, b,)u(R, Fy). (7) 


Since all large integers in turn become central indices of both ¢,(z) and 
F(z), it follows that the values rR include all numbers exceeding some 
definite bound. 

From (5) we get 


A8 


a, (rR)" : n= Rn 


<— 
a,\(rR)’ ~~ vXRY 
P(i+a)(1—R)-* 


HhA8 


f 








py* R’ 
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as R-> 1—0[(2) 225]. Now 


: (aloga+zrlog R) = *+ log R, 
x 


dx 
so that x*R* increases steadily for x > 0 till 

x=£= —a/logR 
and then decreases steadily. Hence, if v = »(R, F,), then |jy—é| < 1. 


Therefore 


l'(1+a)(1— R)-=-1 P(l+a) _ xv { a a — 


~ a~*-!I[P'(1+a)e*.v, (9) 


as v->oo. Since, as has already been observed, all large values are 
assumed by rR, we get (1) from (6), (8), (9). 
Also, from (1), 
. dX , . ; 
lim sup SHU.) <1+lim sup O8¥rf) 
roo logy(r,f) roo logy(r,f) 

and, using (6) and (7), we see that the limit on the right-hand side above 
does not exceegl 


lim su log B, Fa) < limsu logy as 
R-1 Piog u(R, F.) ~~ P vlog yay (v—1) a’ 





which gives (2). 
This completes the proof of Theorems | and 2. 


3. We shall now construct the function g,(z) of Theorem 3. Let 
r, <1) <... <1, be k positive numbers. Define 


r, = r,exp{k-+...4+(n—1)-} 


for k < n < p?, where p, is chosen so that 


rk-lexp4k-3 
k < pi, re. exp gk™ < log p,. 
es el 
Then we get 
“Pi ‘—1 rk-1 } ;-3 
_"n et exp{k-4+...+(p,—1)} < +———*+— exp(s*™) 
TyeM yp, VT yroP gun 1-0F p14 
< log p,. 


3695.2.7 N 
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When p, < n < pi, 


nr Pr n~py 
Bf a 


yf] Typ, Tpira00Fn 


geek O—,  Remeeess >» B*) 
= exp| pe (pt) "* @—y 


one a 
1 


> exp(—pj'), 


and, when p} < n < py, 





n p 
_Tp - |B a Tn + uel p, 
1y-Ty/ yb p, eS” 
», —n—1 », —n—2 ] 
= exp| —|2 42 ae fee tH | 
\(~i—1)® © (p12) (n+-1)°} 


pp) (a2 1") 
> oP EI | 





> exp(—py?). 


Thus, if p} < n < pi, 


-n rPi 
-™\_ > exp(—p;!)—*-. 

pe pelle, 
Now define r,.,, = 7,;+ 1 and proceed as before with k replaced by p}+1. 


Then we obtain an integer p, > p, such that 


rP? 
—"_ < log p, 
po 
and, when p} < n < pi, 
r” Pr 
at ie exp(—pz*)—“-. 
Fane ee me 


Proceeding step by step in this way we construct a sequence of positive 
numbers 1;,.-., 7,;--- Such that r, increases steadily to infinity with n. 
For an infinite sequence of integers 7y,..., Py;-+- 


™ 
inns << SUM 10 
anol ZPn (10) 


Pr 
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and, for pt < m < p?, 
r™ Pn 
4 Re —_ exp( —p;,*) _.. (1 1) 
Vy-+Tm Ty---Tp, 
, ~. m%zn 
We now put g,(2) = > Pash. ION 
Ty 


Since r,, increases steadily to infinity with n, it follows that g,(z) is an 
integral function satisfying the conditions of Theorem 1. From the proof 
of Theorem | we have v(r,, R,,g,) = P,, where 


x a 
ee ed eee: 12 
one — 59) << 9-55) 08) 


From (11), 
i pi 


ms ™ ~ exp(—pzt)rP= 
1dtp Be) > > a > SS ress mR. (13) 





n a 


S ae ee % 
Also 
PS pt 20 
> m*R™ > {1—e(p,)}T(l+a)(1— R) ~a—1__/ >+> )m=Rr, (14) 
Ps _ & 


where ¢(p,,) > 0 as p, > 00. For the sums on the right-hand side of 


(14) we get vt pa? 
= mee a SR 
5 ‘ i—H,, 

ow anit (15) 


as p,, > 00 from (12), and 


2 1c 
> m x Rm =— . (m x mR ym 
-_ n y= n 


P, Py 


where, from (12), 


5) , 
R,, : exp(— il ad < exp( — = } (16) 





z Pyotl” PR 2ePn 
provided that p,, > (2e—5)-!. Therefore, if p,, > (2e—5)-, then 
J RiP, 
> ate 2" 
pt ? n ~ 1— R;, 
< (1-+-0(1)}2a-%ep,exp( — 2), (17) 


from (16). Since, from (12), 
(1=R,)-2- ~ a-*T (1-0) ht 
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as p, > 0, it follows that 


Pra 
Gall'p, Ry) > {1—0(1)}a-*T(1+-a)pht9— "Ps, (18) 
Pyeng, 
from (13), (15), (17). Now 
“ rp R,, a 
rp, By» Gq) = Pate Sal” (19) 
ee 
and, since R?» ~ e-“ as p,, > 0, we get, from (18) and (19), 
g Ar) = {1 =—— ( ! )} Se “EX 1 +a)e%(r, Ja)u(T, Ix) 
as 7 —> 00 through the sequence r,, R,. This proves (3). 
From (19) we have, using (10) and (12), 
Hp, RnsJa) < piexp( — 5s; \ou Pe 
; log p,, ~! 
and so lim sup ~ (20) 





Tp, ltn>2 log u(r, R,.§ qa) a 
Since p, = v(r,, R,,g,), we get (4) from (19), (20) and Theorem 2. 


This completes the proof of Theorem 3. 


4. I shall now show how Theorem 2 can be used to generalize a result 
given by Shah and Singh (1). These authors have proved the following 
theorem. 

Let (x) satisfy the following conditions: 

(i) O(a) is positive and non-decreasing for x > x, and tends to infinity 
with x, 

rd 
‘i at 
ll) I(x) = wears 
( ( 0 
tends to infinity with x, 
x6’ (x) 
Oz) 
for x >x,. Let N be an integer such that I(N) > 1. Then 


(iii) 


F(z) = s {z/I(n)}" 
< 


is an integral function of infinite order such that 


f 7 
rx es F) 





If 6(x) also satisfies the condition 


(iv) I(n®)—I(n) > 





O(x+1) 
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where p is some (fixed) integer for all large n, then 


lim OB Mr, PyOilog M(r, FY} _ 9 (22) 
rw v(r, I ) 





Shah and Singh showed that (22) follows from (21) and (i)-(iii) if 
log M(r, F) ~ logu(r, F). To prove this last result they made use of (iv). 
However, I shall show that log M(r, F) ~ log p(r, F) follows from (i)—(iii) 
alone. 

Let =r 
x(x) =A log( : + (x+ l)log J(x+1)—z2z log I(x). 


Then, using (ii), we have 


A l 
x(x) = - 4 —~ + log I(x+-1)— 


a(x+1) Oa+l)I(a+1) , — — log I(x). 


l 
A(x) 1 (2) 


From the mean-value theorem we get 


] 
log J(x+-1)—log I Sp. onoenenerintemnienmenn ee 
S5e-be (x+-1)0(x+1)L(x+1) 
and, using (iii), 


l l z c 
Oax+ljl(x+1) O(x)I(x) ~ - gaTe- (2) 





] 1 l 
« 1 - —~-— —-—— = = ——j, 
wast (x+1)O(a+1)I(a+1) x(x) (x) - o(5) 
From these results we find that 
S(t) > A F 1\, Il-—e 
saudi ~ a(a+1) x?) ' x0(x)I(x) 


Condition (iii) gives 6(x) < {1+-0(1)}a* and, since J(x) = o(log~), it 
follows that 





> 
= 


= BD me 
ie) > — 4 of) 4 U9, 


a(z+l1) al log x 


which is ultimately greater than zero for any A. Hence x,(x) is ultimately 
a steadily increasing function of x, and so 


n+) {I(n+1)}"*1 
{I(n)}" 


exp{x,(%)} = (' ~ 
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is ultimately a steadily increasing function of n for any A, however large. 
From Theorem 2 we get 


log M(r, F) ~ log p(r, F). 


This shows that condition (iv) is superfluous. 
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A CLASS OF MAJORITY GAMES 
By J. R. ISBELL (Washington, D.C.) 
[Received 16 December 1955] 


Tuts paper is concerned with the description of a class of combinatorial 
objects, the weighted majority games of von Neumann and Morgenstern 
(2). I find a rather full description of a certain subclass containing 2"~* 
n-person games for each n not less than 4. This class is probably a 
negligibly small part of the whole; but the corollary that there are at 
least 2"-* different n-person weighted majority games is better than 
any previously known lower bound. (For six players there are exactly 
14 such games, and for seven at least 110.) Beyond our little subclass 
we have a few remarks; but it is not even known whether in general there 
is a unique natural choice of weights. The homogeneous weights are unique, 
when they exist (1); and I show that such weights cannot increase faster 
than the Fibonacci numbers. 

A (strong, weighted) majority game is a game definable by weights w; 
assigned to the players; the set S of players wins if and only if its weight 
w(S) = > w; is greater than the weight of its complement N—S. 

icS 


The essential fea consists in the specification of winning sets (and their 
complements, the losing sets). There can be no ties. Thus we deal with 
n-tuples of non-negative real numbers, organized into equivalence classes 
which are essentially open convex cones. I shall suppose that no dummy 
players are present, so that each player i is in at least one minimal 
winning set. 

Given an abstract game, i.e. a specification of winning sets, one can 
determine whether it is a majority game by constructive methods: it is 
a question of solvability of a system of strict linear inequalities. How- 
ever, no simple constructive characterization is known. Thus we face 
the questions, ‘What pairs of n-tuples of weights are equivalent?’ and 
‘What properties have the equivalence classes ?’ Since I have no answers, 
I shall minimize the discussion of these questions. But obviously open 
cones contain lattice points, i.e. there are integer weights for any 
majority game; and minimal integer weights are easily found. In all 
known examples there exist minimum integer weights, i.e. the cones in 
question have the very unusual property of containing minimum lattice 
points. Any substantial contribution to the theory of weighted majority 
games ought to determine whether this property is accidental or general. 
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For each player i in a minimal winning set S, there is a minimal 
winning set 7’ whose intersection with S consists just of 7. This is well 
known and follows from the fact that (N —S) U {i} wins but V—S does 
not. Several corollaries follow directly. 

Let us write S’ for N—S. Note the identity 


w(SO T) = w(S8'N T’)+4[w(S)—w(S8’)+ w(T)—w(T")). (1) 
When S and 7 are winning sets with intersection {i}, this means that 
w; > w(S’N T"). 
[1] Each pair of players lies in a minimal winning set. 

Proof. Let the players bei and j (w; > w;). Let S and T be winning sets 
meeting in {j}; since w; > w(S’N 7”), i is not in 8S’ N 7” and hence is in 
S or in T. 

[2] Lf each player is the intersection of two sets S such that 
w(S) = w(S’)+1, 
then the weights w; are minimum integer weights. 

Proof. Let w; be the smallest non-integer weight. From (1) there is a 
set S’NT” of weight w;—1, which is absurd. Then the weights are 
integers; suppose that they are not minimum. Let (v,) be integer weights 


and w,; the smallest weight such that v; << w,;. From (1) there is a set 
S’'N T’ such that w(S’N 7”) = w;—1. By choice of j, 


v; J u,;—l, (SN T") > w(S’NT"). 
Substitution gives v, < o(S'NT"), 
a contradiction. 

The weights are said to be homogeneous when w(S)—w(S’) is the same 
for all minimal winning sets S (2); letting this common difference be 
unity, we see that 
[3] Homogeneous weights are minimum integer weights. 

[4] Homogeneous weights are unique. 

Corollary 4 has previously been established by an argument on the 
rank of the incidence matrix of the game (1). 

From the hypothesis of Corollary 2 we have deduced that for each i 
there is a set of weight w;—1. A simple induction proves that 


[5] Every integer from 0 to w(N) is w(S) for some S. 


With homogeneous weights we can use Corollary 1 to set aside the 
player j whose weight is the largest of those less than w;: that is, i and j 
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are in a common minimal winning set S and hence there is S’N 7” 
excluding j such that w(S’N 7”) = w,—1. This leads to 


[6] The k-th from least weight of homogeneous weights is at most the 
k-th Fibonacci number. 


Proof. Recall the definition of the Fibonacci numbers: 


ay — as _ 1, a, — A, -11+En-2- 


By induction, the sum of the first » Fibonacci numbers is a,,,,—1. Now 
the least weight w; is 1; and the kth from least exceeds the sum of the 
smallest k—2 by at most 1. Hence the corollary. 

With Corollaries 1, 5, 6 (and 3), one has the beginning of a description 
of those sets of integers which are homogeneous weights. Corollary 6 is 
sharp; letting the first n—2 weights be the corresponding Fibonacci 


numbers, w.=a@ 


Wy-1 = Fn-2 n n-1? 
one has homogeneous weights for each n > 3. 

We turn now to the narrow class C of all majority games with n (> 3) 
non-dummy players and nm minimal winning sets. From (1) we know 
(a) that this is the least possible number of minimal winning sets, 
(6) that each game in C has a main solution, (c) and therefore that each 
game in C has homogeneous weights. Further, we assert the next three 
paragraphs. ~ 

The class C consists of the three-person majority game (weights are 
1, 1, 1) and 2"~* games of n persons for each n > 4. Forn > 4, arranging 


the homogeneous weights w; in non-increasing order, we have 
W, > We = Wy, Wn = W,. 
This leaves n—4 choices, 


WO, > Wz, OF Wy = Wis, 
which may be made freely. The players are then partitioned into h 
equivalence classes £,,...,£,, arranged in order of decreasing w;. One 
minimal winning set S, consists of the union of all Z; such that j = h 
(mod 2). The other minimal winning sets S; may be indexed by the 
players i > 1. For? in &,, S; consists of i together with all Z,, such that 
k <j and k = j+1 (mod 2). 

The transpose of the incidence matrix of players upon minimal winning 
sets for a game @ in C is the corresponding incidence matrix for another 
game G* in C. The three-person and four-person games are self-dual; 


ee Bo a or eC i +f an , ; * 
for 3 <i < n—2, w,; > w;,,, in G if and only if w,,_;., > w,-;42 in G*. 
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The minimal winning set S, corresponds to the player | in G*; the other 
S; correspond to n+2—i. 
For i in E,, w; = 1. For i in E; such that 
1 <j = h+1 (mod 2), 

w; is the sum of the weights of all players in E, (k >j; k =h (mod 2)). 
If 1 <j = h (mod 2), then w, is the sum of the weights of all players in 
E(k > j; k = h+1 (mod 2)) plus 1. Finally, w, is the sum of all w; (i > 4) 
plus 1. 

The reader may verify that we have correctly described 2"~* games in 
C for each n > 4. Thus our assertions become justified when we show 
that there are at most 2"-* games of n persons in C. We do this by an 
induction, reducing n-person games to (n—1)-person games and showing 
that the transformation is at most two-to-one. For the initial cases 
n = 3andn = 4, we refer to (2). We remark that besides the involution 
G <> G* and the reduction just mentioned, there are at least two other 
natural transformations in C. One consists of deleting all players in £,, 
i.e. of weight 1; the remaining weights are of course not homogeneous 
but they define a game in C if there are three players left. Another 
transformation gives an h-person game (if h > 3), by aggregating each 
equivalence class into a player having the total weight of the class. 

Our argument begins with a fifth transformation; but this is only 
a hypothetical construction leading to a proof by contradiction. 

[I] For G in C with n > 4 players, each minimal winning set contains 
none, one, or all of the players of weights 1. 

Proof. Suppose [1] is false in G. Let w be the homogeneous weights for 
G, and define G’ by deleting two players, i, 7, of weight 1. The sum of the 
weights in G’ is odd and hence G’ is a strong majority game with n—2 
players. It contains players of weight 1 (or else [I] holds); and, since [T] 
is false, they are not dummies. Hence G’ has no dummies and therefore 
it has at least n— 2 minimal winning sets (1). Now consider the t (> n—2) 
partitions of the players of G’ into a minimal winning set S and its comple- 
ment 7’. Either w(S) = w(T)+3 or w(S) = w(T)+1. In the first case 
S is a minimal winning set in G containing neither i nor j. In the second 


case US; Jy si az 
Sufi}, SUf{j}, Tu {ij} 


are three different minimal winning sets of G. Since G has only n 
minimal winning sets, the second case cannot arise more than once. But 
from n > 4 and the denial of [I], one easily sees that this is a contradic- 
tion. Therefore [I] holds. 

















A CLASS OF MAJORITY GAMES 
[II] Each player of weight 1 is in just two minimal winning sets. 

Proof. We may begin by noting that [IT] holds in the unique example 
forn < 4. Now assume n > 4. Delete one player i of weight 1. There 
remains a set of players with integral weights whose sum is even; but 
from [I] we know that, if this set is partitioned into two complementary 
sets of equal weight, then all players of weight 1 are in one element of 
the partition. Thus by replacing each weight 1 by 1+, for sufficiently 
small positive e, we obtain a majority game R(G@) with n—1 players who 
are obviously not dummies. As before, R(G@) has at least n—1 minimal 
winning sets. Now the minimal winning sets of G which contain i occur 
in pairs, S, S’U {i}. For each such pair, exactly one of S and S’ wins 
in R(G@). All other minimal winning sets of R(G@) win in G, and therefore, 
(since G has only one more minimal winning set than R(@)), 7 is in only 
two minimal winning sets, as asserted. 

Now [II] is proved. Note that R(G@) must have just n—1 minimal 
winning sets; therefore R(G@) has homogeneous weights, and, if we 
substitute these for the weights just constructed, the reduction can be 
iterated. 

How is G@ obtained from R(G)? At the risk of boring the reader, I 
emphasize that the weights given for R(G) in my construction are not 
given with the game. What we do know is that some minimal winning 
set S of R(G@) does not win in G; an extra player 7 is added, and SU {i} 
becomes a minimal winning set. One property of S which we may deduce 
from the weights for G is that, if any player j not in S is added to S, then 
the result SU{j} contains at least one minimal winning set different 
from S. Now R(@) is a homogeneous weighted majority game; and this 
property for its minimal winning sets is easily seen to be equivalent to 
containing at least one player of weight 1. Since R(@) is in C, we know 
from [T], [II], and the substitutability of players having the same weight 
that there are only two such sets inequivalent under motions of the 
game. Therefore, there are at most two n-person games G in C having 
the same R(G); and there are at most 2"-* n-person games in C, for 
n > 4. 
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‘ALMOST CONVERGENCE’ AND UNIFORMLY 
DISTRIBUTED SEQUENCES 
By G. M. PETERSEN (Swansea) 


[Received 31 January 1956] 


1. I sHatu describe a method of summation a by a rectangular array of 
linear transformations of a sequence {s;}; 


i ne -. 
(2) (2) (2) 
i 
(m) (m) (m) 
~~. . . 4 


© 
where t” = ¥ al’s,;. I shall say that a sequence {s;} is «-swmmable 
n hn 4 © Cu 


if the transformations ¢!”") converge uniformly in (m) to a limit. One 
method of summation of this type is the method of ‘almost convergence’ 
[Lorentz (1)]. If we let 


(m) __ —-l/(oe ; 
t =m (8, +841 + +->+8nim) 


we have ‘almost convergence’ described by the uniform convergence of 
the rows of the corresponding array. It is also evident that any matrix 


method can be described as a method of this type. If A = (a,,;) and 


mi) 
ms 


tn = > Ang S;, We need only take t!” = t,, for all n. 


m — mi § 
i=1 


Associated with each method of summation « is the family Y of all 
possible matrices that can be developed by selecting the first row of 
the matrix from the first row of the array, the second from the second 
row of the array and so on. This leads to the first theorem. 


THEOREM 1. The method « sums exactly those sequences summed by 
every member of U. 

Proof. If a sequence is not limited by some member of UY, then it is 
evident that the rows cannot converge uniformly. Hence « is certainly 
contained by all the members of Y. 

On the other hand, suppose that a sequence {s,} is limited by all the 
members of YU to s but not by a. Since it is not a-limitable, there is an 
infinite set {m,} and a corresponding {n(m,)} such that 


(my) 
bximy) —? >€ 
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for some «. This means that the matrix which has an infinite number 


of rows selected from fe6mu)) 
n(™my)I? 


and an infinite number of rows from some matrix that limits {s,} will 
be a member of Y& that does not limit {s,}. This contradiction proves 
the assertion. 

We see incidentally that « is regular if and only if all the members 
of U are regular. The method can be strongly regular in the sense of 
Lorentz (1) if and only if all the members of & are strongly regular. 


2. We suppose that 0 < s, < 1 for every n, and denote the interval 
0<a<2<b<1byTl. We denote by /(x) the characteristic function 
of J which is 1 in J and 0 elsewhere. 

A sequence {s,} is then said to be uniformly distributed if 


or ae 
Serr 





for every J. 
A well-known theorem due te Wey] (2) states that, if {s,} is uniformly 
distributed, then 


lim —— a2 f(s.) = | f(x) dx 


n—-x n+1 


for every Riemann-integrable f(z). 
In this note, I shall mean by a well-distributed sequence {s,} one for 


lim — > I(8,) 


p-o Pp fad 


which 


holds uniformly in n. In other words, {s,} is well-distributed if {J(s,)} is 
almost convergent to b—a for every J. Since any almost convergent 
sequence is summable (C,1), a well-distributed sequence is necessarily 
uniformly distributed. However, the converse is not true. 

Let us suppose that a sequence {s,,} is uniformly distributed on the 
interval (0,1). From this sequence we construct a new sequence by 
letting s’, = 1 when n is of the form v*+ 1.,,..., 3+ for all vy and s,, = s,, 
otherwise. This new sequence {s‘,} is also uniformly distributed, for, if 


n 
an interval excludes or includes the point 1, the value of > J(s,) can 
K=1 


be changed only by o(n). On the other hand, for any interval (a, 1) 


viv 
> I(s,) = 1, 
y 


<i 
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and {s),} is not well-distributed. Hence, we see that not all uniformly 
distributed sequences are well-distributed. 


THEOREM 2. If {s,} is well-distributed, then 
1 ; 
lim — &) = d 
ie saa jf 


holds uniformly in n for every Riemann-integrable f(x). 


Proof. The proof of this theorem closely follows that of the corre- 
sponding theorem for uniformly distributed sequences. It is true for 
any finite step-function; two such step functions are selected so that 
fy <f <fy and 
0<|f pdx —[f,dx<e. 

If we form any matrix A = (a,,,) by extracting for our nth row one 
of the transforms in the nth row of our array for ‘almost convergence’, 
then 

co ox - 

lim inf > me f(sx) > lim 2, anetilse) = [ fale) dz = | f(x) dx —e. 

<= mo k= Sf - 


m—>o 


lim sup s And (8%) < [ f(a) dx +e. 
Ko : 


m—>2 


Since « is arbitrary, we have 


lim > Ani f (se) os [fe dx. 


m2 k= 


Since this is true for any of the matrices belonging to ‘almost con- 
vergence’, by Theorem 1, {f(s;,)} must be almost convergent to. { f(a) dz. 
It is now a corollary of a theorem of Lorentz (1) that, if {s,} is well- 


distributed, then in , 
lim a Onkf (8x) ey | f(x) dx 


mo 
for every Riemann-integrable f and for every regular matrix for which 
a 
lim > Amk —4m,k+1! —_ 0. 
m—-x2 k=0 


3. I proceed to develop my criterion for well-distributed sequences. 


ad f(x) = eri — efx), say, 


where k is a positive integer, then { f(x) dx = 0. It follows that {e(ks, )} 
is almost convergent to zero if {s,,} is well-distributed, i.e. {7(s,,)} is almost 
convergent to zero, where 7' is any trigonometric polynomial without 
constant term. 
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THEOREM 3. The sequence {8,} is well-distributed if {e(ks,,)} is almost 
convergent to zero for every positive integral k. 


Proof. Since ‘almost convergence’ is linear, we see that, if 
N 
r(2) = Yagt ¥ {ape(kx)-+b,e(—kz)}, (1) 


then {r(s,,)} is almost convergent to 4a). Now, if J(x) is the characteristic 
function of any given interval (a, 6), we can find functions 7,(x), 7,(x), 
both of the form (1), such that 
1 
7,(x) > I(x) > 7,(2), | [7,(x)—7,(x)] da < e. 

0 
We can thus show by an argument similar to the proof of Theorem 2 
that {/(s,,)} is almost convergent to b—a. The proof can now be carried 
out in the same way as the proof of Theorem 2. 

If s,, = n€—|né| = {n€}, where é is irrational, then it is well known 
that {ng} is uniformly distributed. I shall show that {n€} is well-distri- 
buted. Indeed 
l — e2kip+irgi 


a af k ow 2kuréi 2kn7géi 
8.) = pepets — ete me 
7 $,,) en ¢ 1 —e2k7éi 

= O(1), fork = 1, 2, 3,... and for all x and p. 


Hence {né} is well-distributed. 
I am grateful to the referee for his comment concerning the definitions, 
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ON MEANS OF ENTIRE FUNCTIONS 
By Q. I. RAHMAN (Aligarh) 


[Received 21 February 1956] 


1. Let f(z) be an entire function of order p and lower order A. Then, 
5 and «x being any positive numbers, let 


Pr. 
pg(r) = = f(ret®) > dé, 


ms ,(r) = — ae i f(xet®) |Sxx dxdé. 
0 


Vlog r 
Ls, = lim apf Hal") | 
. ro (ms_(r)} 


Vlog r 
ls, = lim ing |_H8(”) I : 
rw lms «A Ar) 
It is known [(2), problem 66; (3), 96-99] that for every entire function 


L,, = e and 1,, = e*. I prove here the theorem: 


THEOREM. For every entire function Ls, = e?, ls, = &. 


2. Lemmas 
Lemma 1. Jf f(z) is regular in |z| < R, and if 





z= ret? 0<r<R, 8 > 0, 
| 1 (R?—r?)| f( Re’) | 
th yr< = | 18. 
i f@)" < 2r | R?—2Rrcos(—¢)+r°- (1) 
0 


Proof. Case (i). If f(z) has no zeros in |z! < R, then any branch of 
{f(z)\® is regular in |z| < R, and we have, by Poisson’s formula, 


{f(z)}> = = [ (R?—r?){ f( Re'#)}? 
“? “29 | R?—2Rrcos(@—¢)+r? 
3 





from which (1) follows. 

Case (ii). Let f(z) have zeros a,, a,..., «, in |z| < R but no zero on 
'z| = R. Then, using the function which represents |z| < R on itself 
and carries z = 0 to a; (i = 1, 2,..., p), we can construct a function ¢(z) 
having the same zeros as f(z) in |z| < R such that |¢(z)| < lin |z! << R 


Quart. J. Math. Oxford (2), 7 (1956), 192-5. 
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and |¢(z)| = 1 on |z| = R. Applying Case (i) to {f(z)/¢(z)}° we get the 
result in this case. 

Case (iii). If f(z) has a zero on |z| = R, we can establish the result 
for |z R-+. and let «> 0. 

So the result is true in all cases. 


LEMMA 2. Let 





A lim sup OB/OB Hal") B= lim inf 08108 Ha") 
r>x log r r—>2 logr 
x = limsup loglog ms,.(7) | B = lim ing OBB Mo(7)_ 
ro logr rao log r 
Then A=a= Pp, B= B = A. 
Proof. If M(r) denotes the maximum modulus of f(z) for |z| = r, then 
27 
pg(r) = pA f(ret®) 5 dd < {M(r)}. (2) 
aT 
0 
From (1) we have 
R+r 
18) )5 — pe Th ltt R 7 
f(ren)" Se bal ®) 
Taking r = }R, we get 
: us(R) > HM (ER). (3) 


From (2) and (3) it follows that A = p and B = A. 
Since ps(2) is (1) an increasing function of 2, 


r 














2 
ms,,(7) = —> [vse dx 
he 
0 
2us(r) 
9 ° ‘ 
_ 23(r) __ 2palr 
ce eee 
0 
and we get a <A, B< B. (4) 
Further er 
9 
ms ,(2r) = (Ore | ps(x)aX dx 
0 
2r 
2 
> [ ps (ax)aX dar 


(2ry* J 
r 


2pus(r) (2r)*+t—rxtt 
a (2r)*+1 «+1 


3695.2.7 oO 





> 
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and we get a >A, B> B. (5) 
Thus, from (4) and (5), «a = A and B = B. 
Lemma 3. [fd and « are any positive numbers, r**y3(r) is a convex 
function of r<*Ims Ar). 
Proof. We have 
d(ra*usr)) _ (aldr)(re*uslr)) 
d(r**4ms(r))  (d/dr)(r***ms,_,(r)) 





__{k + 1 )r*ps(r)+r***u5(r) 





(d/dr)((1/x7) { { \f(xet®) Sx dard®) 
06 


(«+ 1 )r*us(r)+r*+ty5(r) 





27 
2r*(1/27) { |f(xei®)|> do 





0 
_ (e+ )rXpa(r) +18 ty5(r) 
2r*yus(r) 
1 H3(")) 
= —{(«+1)+r—}, 
z| ua(r)) 


which increases with r since, by (1), log us(r) is a convex function of log r, 
and the lemma is proved. 


3. Proof of the theorem 

If Ls, < 0 then, for a positive « and a suitable constant a, 
r 
| r*5() 


log{r**+4ms (r)} = | —— 
_ nee" ats (x) 





| 


f pa(x) dx 


ms A(x) x 


r 


O()+ | (Ls,,-+e)loex 
*, fhe? 

< o(1) + aatel?8" =... 

log( Ls .+€) 


/\ 


(6) 


From the fact that logys(r) is a convex function of logr it follows 
very easily that 


lim log ms,,.(r) 
r+o logr 


= ©. 
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Hence, from (6), we have L;,. >, 
which obviously holds when L;,, = 00. 


It follows from Lemma 3 that y5(x)/ms (xz) is an increasing function 
of x and therefore, for 0 << Ls, < a, 


yy 


_ @*yg(Z) _ 


log{(2r)*+4ms .(2r)} > j 
a a**lms (x) 


2r 
Balt) [ dx 


- so 
ms <r) J x 
r 


_- halt) lo 
ms (1) 


> (Lz,,,—«)'°€" log 2, 





for a sequence of values of r tending to infinity. Consequently 
D5. S 
which holds when Ls, = 0. If Ls, = 00, the above argument gives 
e? = 00. This proves that Ls, = &. 
The proof that /;,, = e is similar and is omitted. 
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THE MINIMUM MODULUS OF FUNCTIONS 
REGULAR AND OF FINITE ORDER IN THE 
UNIT CIRCLE 


By C. N. LINDEN (Cambridge) 
[Received 19 March 1956] 


1. | SHALL examine the association which exists between the minimum 
and maximum moduli of certain functions regular in the unit circle. 
The general results which I obtain here are most significant when applied 
to a function g(z) of finite order 





; loglog M(r,q) 
= lims 3 a 
B im sup =i 


where M(r,g) = max'g(z)|, the minimum modulus function 
=r 
p(r,g) = min'g(z)| 
2|=r 

being compared with M(r, g) for a sequence of values r which increase to 1. 

The corresponding problem in the theory of integral functions has been 
considered in various ways. Littlewood (4) proved that, if G(z) is a non- 
constant integral function of finite order p, then 


. lo og p(T, G) 
parse fy ae: > —C(p), (1.1) 


where C(p) depends only on p. For large values of p, Hayman (2) has 


shown that C(p) < 2-19log p, 


and that, if G(z) has infinite order, (1.1) must be superseded in general 


by 
> —2-19. 





lim s sup log u(r, @) . 

r+o log M(r, G)logloglog M(r, G)~ 

The method used by Hayman to prove these results is extended here to 

apply to functions regular for |z| < 1. The modifications which I make 

in the theory take the form of a number of preliminary lemmas which 
generalize this earlier work. The main theorem which I prove is 


THEOREM |. Let g(z) be a function regular and of finite order B for \z| <1. 
Then, if 8B > |, there is a corresponding constant K(B) such that 





—K(8) < limsup log (7, 9) < 0. (1.2) 


I 


r+1 log M(r,g)loglog M(r,g) 
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The symbols K, A, B, ete., which arise in the sequel represent finite 
positive constants, not always taking the same value at each occurrence 
and sometimes dependent on certain parameters which will usually be 
indicated. 

Theorem | differs from the result quoted above for integral functions 
in the occurrence of the term loglog M(r,g) in (1.2). The necessity to 
include such a term in a general result of this type is shown later when, 
for each finite value 8 greater than 1, we construct a function g(z) regular 
and of order 8 in |z| < 1 for which the limit occurring in (1.2) is non-zero. 
[ remark also that the inequalities (1.2) are not true in general for 
functions which are regular in |z| < 1 but of order not exceeding the 
value 1. This type of function is not dealt with in any great detail here 
except to consider briefly how the behaviour in particular cases differs 
from that of general functions of greater finite order. 


2. The earlier lemmas of this paper are proved with respect to a 
function f(z) regular for |z} <1 and an associated function F(w), 
extensions being made which apply to functions g(z) regular for |z| < 1. 
With Hayman I write the function u(z) = log|f(z)! as 


u(z) = u(z,h)—N(z,h,f). (2.1) 
‘ h 

where N(z,h,f) = n(z,t,f) t-' dt, (2.2) 
0 


and n(z,t,f) is the number of zeros of the function f(£) in the circle 
z—&| <t. Thus u(z) is expressed as the difference of two components 
by considering those zeros of f(€) which are situated in a certain neigh- 
bourhood of the point z. I shall examine each of these terms separately, 
employing the following two lemmas proved by Hayman [(2) 476, 483]. 


Lemma A. Suppose that F(w) is regular for |w| < 1, F(0) = 1 and 
0<p<1,h = &(1—p), where 0 <8 < }. Then, if U(w) = log|F(w)!, 


we have 


27 
U(p,h) ee | U +(e) dy, (2.3) 
0 
27 
n(p,h, F) <i +2 t | U+(el4) dy, (2.4) 
l—p 22 , 


where U+ = max(U,9). 
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Lema B. Suppose that 0 <r < 1,h > Oand that for |z| = r we have 
n(z,h) < m. Then there is a set of values o for which r <a <r+th 
having measure at least th and such that 


N(z, fh) < nolog >, 


for |z| =o. 

In Theorem 2 we deduce inequalities for the functions u(z,h) and 
n(z,h,f). Thence, if g(z) is regular for |z| < 1, we obtain a corresponding 
bound to u(r,g) in Theorem 3 in terms of an integral which involves 
M(r,g). The result is brought into the particular form (1.2) by intro- 
ducing the concept of ‘proximate orders’, which is familiar in the general 
theory of integral functions. 


3. Functions regular for |z| < 1 

By taking the integral of U’+(w) over the whole of the circle |w| = 1 
in Lemma A we may, especially if U +(w) is large on an appreciable portion 
of that circle, fail to obtain bounds to |U(p,h)| and n(p,h, F) which lead 
to results of the best possible order. Roughly speaking it is only that 
contribution to the integral arising from an appropriate neighbourhood 
of ys = 0 which is necessary to fix the bounds to |U(p,h)| and n(p,h, F). 
Consequently we make a transformation which effectively reduces the 
range of integration of the integral which occurs in (2.3) and (2.4), thus 
sharpening the applications of these inequalities. Certain properties of 
this transformation are given by Lemmas 1, 2, 3. 

Lemma 1. The transformation 

1—z = (1l—w)*, 4$<a<l, (3.1) 
maps the circle \w| < 1 onto a region G which, except for the point z = 1, 
is entirely contained in |\z| < 1. 

Hayman has previously considered the particular mapping given by 
(3.1) with a = 3. We examine that branch of the right-hand side of 
(3.1) which is positive when | —w is positive. Thus G is simply-connected 
and each of its boundary points corresponds to some point w = e' 
(0 < % < 27) of the circle |w| = 1. For such a point z we have by 
(3.1) 


z= 1—(1—el+)* = 1 —eliald—2(2 sin 4h): 
that is z|? = 1—2 cos $a(y—7)(2 sin $y)*+ (2 sin 4y)**. (3.2) 
In order to find an upper bound to |z! we examine 


I(b) = 2cos $a(—7)— (2 sin y)*. (3.3) 
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The function J(y%) is continuous in the interval (0,27), taking at each 
of its end points the value 


I(0) = I(27) = 2cos $a = 2s8in }(l—a)mr > 2(1—a), (3.4) 
since } << «<1. At the turning points of J(%) we have 


dl 
dds 


so that for these values of % in the interval (0, 27), 
2sin $¥% sin $a(y—z7) 


—a sin da(y—7)—a cos $4(2 sin $y4)*-! = 0, 


T(x) 2 cos $a(ys—2) + 








cos $y 
, cos {y—(p—m)a} sin }(1—a)(wm—p) _ 
- cos $y ~ " gind(a#—p) dialing 


Thus from (3.1), (3.2), (3.3), (3.4) we deduce 

z|? < 1—2(1—a)(2sin })* < 1, 
for w = e¥ (0 <& < 22) with equality holding only if 4 = 0. This 
completes the proof of Lemma 1. 

4. We now consider a function f(z) regular for |z| < 1, associating 
with it a function F(w) to which we can apply Lemma A. Lemmas 2 
and 3 contain various relations which exist between the particular 
functions f(Z) and F(w). 

LEMMA 2. Suppose that f(z) is regular for \z| < 1 and that 

F(w) = f{z(w)} = f{1—(l—w)3,_ Sa <1, (4.1) 
u(z) = log|f(z)|, U(w) = log|F(w)}. 
Then to each positive value € there corresponds a value ry = ro(a,€) < 1 
such that 


U+(e¥) db < sd a (cos $ma)-¥* | M(t, u+)(1—t)¥*-1 dt + M(ro, u*). 


2 . arr 


. 


0 ° (4.2) 
Each boundary point w = e of the circle |w| = 1 corresponds by 

(3.1) to a point z where 
z = rel? — a(y) = 1—(1—el#)*> (0 <p < 2n), (4.3) 


and we can write 


27 


[ U+(e%) dp = [ w*[a(p)]+w*fa(27—p)] dy. (4.4) 
0 


0 


oes 
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Considering this latter integral we note that, when 0 < » < ¢ < 72, 
\a(h)| = |a(27—)| < rq = r9(a,n) < 1, 
from which it follows that 
ut[a(p)]+u*[a(2n—yp)] < 2M (ro, w*). 

Hence 
-¢ M(ro, ut 
3y | weladh)-+ tala —yy ay < iro. aa bad < M(rp, u*), 

7 (4.5) 
and it remains to consider the contribution to the integral (4.4) which 


arises from the interval (0, »). 
By (4.3) and a rearrangement of (3.2) we have 


]—r? = (2sin $)*%{2 cos $a(—7)—(2 sin $4)%}, 


and by differentiating we have, as > 0, 





SF — (a cos $7a)—1(2 sin $4)!-. 
. 
Therefore, if 7 = (a, €) is a sufficiently small positive value, 
| 
tI <- € (cos $ara)-¥*(1 —r)¥a-1, 
dr| ~ a 


for 0 < % < ». Since 


u*[a(ys)]+ut[a(2x7—)] < 2M(r, ur), 
it follows that 


7 
ay | wily +w'ta(2—yy ap 
0 


l+e 


ar 


(cos $2a)-"% M(r, ut+)(1—r)¥=—1 dr 
Ja(n)| 
1 
< l+e (cos Jor) | Mit, ut)(1—t)¥a-2 dt. 


aT 


< 





0 
Together with (4.4) and (4.5) this latter gives (4.2). 


5. Although the function F(w) is not regular for |w| < 1, its only 
singularity in that region is at the point w = 1, where it is bounded. 
Hence we can apply Lemmas A and 2 to obtain an upper bound to 
n(p,h, F) and |U(p,h)| in terms of M(t,u+). Corresponding bounds to 
n(r,h, f) and u(r, h) can be deduced by application of the following lemma. 
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Lemma 3. Let the functions f(z) and F(w) be defined as in Lemma 2 and 
suppose that0 <r <1, 





1—p = (1—r)¥e, (5.1) 

O<h<}I—r), hky=FA(—r™4, hy = BA(I—r) 1, 
Then n(p,hy, F) < n(r,h,f) < n(p, he, F), (5.2) 
U(p,hy) < u(r, h) < U(p, hg). (5.3) 


The inequalities (5.2) are a consequence of the next lemma which 
will be proved later in this section. I use the notation of Lemma 3 in 
the statement. 


LEMMA 4, The transformation (3.1) maps each circle 
lz—r| <h < }(1—1r) 
onto a region G(r,h) of the w-plane which is contained in the circle 


lw—p| < hy 
and contains |w—p| < hy. 


We can now prove Lemma 3 in the same way as Hayman proved a 
corresponding result [(2) 480, Lemma 3 (ii)]. The number n(r,h,f) of 
zeros of f(z) in |z—r| < h is equal to the nuraber of zeros of F(w) in 
G(r,h). Sincethe region G(r, h) is contained in |w—p| < h,, by Lemma 4 


we have - 
n(r,h,f) < n(p, hy, f). 
The remaining inequality of (5.2) is true similarly. 
Now by (2.2) we have 
h 1 


N(z,h) = ( n(z,t) tdt = [ n(z, wh) pd, 
0 


and, by (5.2), 

n(p,uh,, F) < n(r,ph,f) < n(p, pho, FP), 
for 0 <p <1,h < }(1—r). Dividing by » and integrating with respect 
to this parameter over the interval (0,1) we obtain 

N(p,hy, F) < N(r, hy f) < N(p, hg, F). 
The inequalities (5.3) now follow from (2.1) since u(r) = U(p) by 
definition. 


In order to complete the proof of Lemma 3 we need to verify Lemma 4. 
By (3.1) and (5.1) we have 


w—p = (l—r)¥*—(1—z)¥*, (5.4) 
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When |z—r| = h < }(1—r), we have 


|\l—z! < 1—r+|r—z] < }(1—7), 


ldw l , 
IS - =—|i—s/""* < 3(1—r)"o-1, 
il 


since } < a < 1. Hence (5.4) gives 


dw! 


lw—p| <|z—r| max < 3h(1—r)"*-! = hg, 


je—-ri<d(1—r) 








~ 
~ 


which shows that G(r,h) is contained in |w—p| < hy. 
To show that G(r,h) contains |w—p| < h,, we prove that the inverse 


transformation 
' . z= 1—(l—w)* 


maps the circle |w—p| < h, onto a region included in |z—r| < h. By 


(3.1) and (5.1) we have 


z—r = (l—p)*—(1—w)*. 


Now h, = $h(1—r)"2-! < 3(1—p), 
so that, when |w—p| = h,, we have 
l—w| > 1—p—|p—w| > }(1—p), 
| dz 


aj = a|l—w|*-! < }(1—p)*-. 
dw} 
Hence, for these values of w, 


| 


~ | 


z—r| < |w—p max j|——| < $h,(1—p)*"' = th < h, 
w—p|<a(l p)|dw} 
and Lemma 4 is proved. 
6. Lemmas 2 and 3 provide with Lemma A the means of proving the 


_— 


main result concerning functions regular for |z| < 1. Itis stated thus: 


THEOREM 2. Let f(z) be a function regular for \z| < 1, f(0) = 1, and 
suppose that} <a <1. Then to each positive value « there corresponds 
a value ry = ro(a,€e) < 1 such that for \z| =r and h = 8(1—r), where 
0 <6 < 4, we have 


1 
u(z,h) > — Aa, 8)(1—r)-14 | M(t, u+)(1—#t)"*— dt + M(ro, vw), 


; (6.1) 


n(z,h,f) < (1+68)(1—r)-#* x 
1 


x pas (cos $77a)-"% | M(t, u+)(1—t)¥*-1 dt + M(ro, w (6.2) 
0 


QT 




















ON THE MINIMUM MODULUS 203 


By applying Lemmas A and 2 to the function F(w) defined by (4.1) we 
obtain 


, _ A(a, 84) 
U(p,h,) > - = {i M(t, ut+)(1—t)"*-! dt + M(rp, w)}, 
for 0 < p< 1, h, = 8,(1—p), 0 < 8, < 3. Now, if 
h = 2h,(1—r)!-¥@ = 28,(1—p)(1—r)!-¥@ = 28,(1—r), 
the hypotheses of Lemma 3 are satisfied for 8 = 25, < }, and it follows 
by (5.1) and (5.3) that 
u(r, h) > U(p,h,) 
(a, 5)(1—r) al f M(t, u+)(1—t)¥=— dt +M(ro,w}, 
\3 
for h 5(1—r). 
Further from (2.4), (4.2), (5.1), and (5.2) we can deduce that, for 
h, = 8,(1—p) < 4(1—p), h = th,(1—r)-¥* = 38,(1—r), 
n(r.h,f) < n(p, ho, F) 


1 


. (cos $7a)-"* | M(t, u+)(i—t)'*—! dt + M(ro, u .)\ 


} 


(1+28,)(14 


l—p | an 
0 


(1+ 28,)(1—r)-"* x 
1 
(1+ ’ | 
-£ (cos hora)-"™ | M(t, ut)(1—t)¥*-! dt + M(ro,u*))}, 


| an | 
0 


provided that 5» < 4. If we choose 
6 = }5, < }, 


the inequalities (6.1) and (6.2) follow with z = r. That the inequalities 
are true in general can be proved by a similar consideration of the 
function f(ze*). 


7. Functions regular for |z| < 1 

The theory of the preceding sections can be adapted to apply to 
functions regular for |z| < 1 by considering these functions in the circles 
z| < R< 1. In this way we extend Theorem 2 and effect the proof 


of Theorem | by making use of certain fundamental results in the theory 
of proximate orders. 
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THEOREM 3. Let g(z) be regular for |z| < 1,9(0) = 1, and suppose that 
$<a<l. Then, ifr < R <1, there is a set of values + of measure at 
least 3;R(1—r) in the interval {Rr, R(r+-34(1—r))} and a value 

ro = (a) <1 


such that, when ry <r < 1, 


v(z) = log}g(z)| > - (1—r)-1 Jog 4) x 
l—r 
> R 
x | (eos ba Va | M(t, v+)(R—t)¥*— dt +2M(ro, v* »\, (7.1) 
X7T 
E ' 
on each circle |z T. 
Since the function f(z) = g( Rz) (7.2) 
is regular for z| < 1, we have, by applying Lemma B and Theorem 2 


with h = 8(1—r), « = 8 = }, that in the interval (r, r+ 4h) there is a set 
of values o of measure at least }/ for which 


N(z, 3h, f) < (1—r)-"* log - lee 


l—r 


< {2 (cos $am)-"* | M(t, u*+)(1—t)¥*-1 dt +2M (ro, uw), (7.3) 
|am 7 , | 
0 
for |z| = o. Hence for such values of z we have, by (2.1) and (6.1), 
log! f(z)| = u(z, 4h) —N(z, gh, f) 
A(«) 
. —lja te 
> —(l1—r)-"*log oe 
; 1 
x (2 (cos 2ra)-"% [ M(t, ut+)(1—t)"*— dt 42M (ro, u*)}. 
S | 
0 
From (7.2) it follows that for |z} = Ro = 7 we have 
log|g(z)| = log| f(zR-*)| 
=) 


> —(l—1) "hg — 


1 
x [2 (cos $7ra)- um [ate u)(1—t+)¥*-1 dt +2M(ro,ut)}. (7.4) 
aT 
0 
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But, since u(z) = v( Rz) by definition, it follows that 
M(t,ut+) = M(Rt,v*). 
Hence we have 


1 1 
{ M(t, u+)(1—t)¥o—- dt = [ M(Rt,v+)(1—ty"=— de 
0 0 
R 
= R-W« { M(t, v+)( Rt)" dt. 
6 
The inequality (7.1) now follows from (7.4) since M(ry,u*+) < M(ro, v*). 
The set of values o for which (7.3) holds on |z! = o has measure at 
least 45(1—r) in the interval {r,r+-}8(1—r)}. Because we have chosen 
5 = }, the set of values 7 for which (7.1) holds on |z! = + has measure 
at least ,; R(1—r) in the interval {Rr, R(r+-4,(1—r))}. This completes 
the proof of Theorem 3. 


8. Theorem 3 is a general result which is valid for all functions regular 
in the circle |z| < 1. We now deduce Theorem | as a particular applica- 
tion to functions of finite order greater than 1. It is required to associate 
the integral which occurs in (7.1) more directly with the function M(r, g) 
and we proceed as follows. 

A function. p(x) is said to be a Lindeléf proximate order if it satisfies 
the conditions [(1) 54] 

(i) p(x) is real-valued, continuous, and differentiable in adjacent 
intervals for x > Xo; 

(ii) lim p’(x)rlogxz = 0; 
at 


(iii) lim p(x) = p (0 < p< @). 


Proximate orders provide a means of approximating to the maximum 
modulus of integral functions of finite order, and it has been proved that 
to each such function G(z) there corresponds a proximate order p(x) such 


that lim sup r-*” log m(r, G) = 1. 


In proving this result, Shah (5) has not needed many special properties 
of the function M(r, @). In fact he has based his proof solely on the con- 
tinuity of M(r, @) and the fact that for functions of finite order, 
lim sup loglog M(r, @) 
roo log r 
is finite. Thus his proof covers more than he has actually stated and 
in particular 
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THEOREM 4. Let (x) be a positive continuous function defined for 
x > 2X», such that 


. log ‘¥'(a) 
lim sup — -= p< 0, > 0. 
Ee ~ P log x > p 
Then there is a Lindeléf proximate order p(x) such that 


lim sup 2-P""F(x) = 1. 


mn 
As a corollary we have a corresponding theorem for functions regular 
in the unit circle. 


THEOREM 5. Let g(z) be regular, unbounded, and of finite order for 
z| <1. Then there is a Lindeléf proximate order p(x) = B(x~*) such that 


lim sup(1—r)P@—) log M(r,g) = 1. 
r—1 
For the function 
¥(x) = log M(1—2z-,9) 


satisfies the hypothesis of Theorem 4. Hence there is a Lindel6f proxi- 
mate order p(x) such that 


lim sup xP) log M(1—a-',g) = 1. 


If we write x = (1—r)-!, B(x-") = p(x), we have 


lim sup(1 —r)pa log M(r,g) = 1, 
r—1 


thus verifying Theorem 5. 
A further property of proximate orders which will be required later 


is that, for each positive constant k, 
} 


(kar)ee) ~ hex), i (8.1) 
as x tends to infinity [(1) 55, Lemma 2]. This can be deduced from (i), 
(ii), (iii). 


9. If p(x) is a Lindeléf proximate order associated as in Theorem 5 
with the function g(z) of Theorem | and 


B(t) = p(t-), (9.1) 
then it follows from (iii) that the order of g(z) is by hypothesis 
lim sup 8(1—r) = limf(1—r) = B > 1. (9.2) 
r—1 r—l1 


Before making our application of Theorem 3 we proceed to replace the 
integral which occurs on the right-hand side of (7.1) by a more con- 
venient expression by writing R = r. 
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Theorem 5 shows that to each constant C greater than unity there 
corresponds a value ry such that 


M(r,v+) = log M(r,g) < C(i—r)-8@ (rg <r < 1). (9.3) 


If r, is chosen appropriately, then by (iii), (9.1), and (9.2) the function 
fef. (1) 55, Lemma 1] (1—r)K8+0-Ba-» 


increases with r for r,s < r < 1. We now define 


1 [18+8) (1<B <5), 


a (65<B<o), 
and note that } < a < 1. We have also 
x? < }(3+8) < H1+8) (lL<B <0), (9.4) 
so that, by (9.3), 
i Mit, y? \(ir—t)? x—1 dt <€ [ (1—t)¥=-1-Ba—» dt 
< C(1—r)i4+8)-Ba--) ( (1—t)! x 1(3+8) dg 
To 


< K(B)(1—r)¥!=-Ba-», 


if, for example, we take C = 2. Further 
( M(t, v+)(r—t)¥2-! dt < M(ro, vt) | (r—t)¥*-1 dt < M(ro, v*). 
0 0 
Thus for a suitable value r, we have 
{ M(t, v+)(r—t)¥—1 dt < K(B)(1—r)¥*-80-+ Miro, v+), (9.5) 
0 
een =F < I. 
Now, if R = r, the interval 
{Rr, R(r+4(1—r))} 
of Theorem 3 is contained in (2r—1,r) for sufficiently small values of 
i1—r. Hence the application of this theorem gives that in each interval 
(2r—1,r) there is a value 7 for which 


log u(7,g) > —K(B)log #. {(1—r)-80—) + M(rp, v+)(1—r)-¥}, (9.6) 


by (9.5). Here we assume without loss of generality that g(0) = 1. 
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By Theorem 5 we have that for each constant D less than unity there 
is an infinite sequence of values s increasing to 1 such that 


log M(s,g) > D(1—s)-Fa-), (9.7) 


We may write any typical value s as 2r—1, and have that (9.6) holds 
for a value 7 in the interval (s,}(1+-s)). Since 
l—r < 1—s = 21-7), 
we have log M(r,g) > log M(s,q), 
and it follows from (8.1) and (9.7) that, if 1—s is sufficiently small. 


log M(r,g) > D(1—s)-®!-* > K(gp)(1—r)-f—”, 
loglog M(r,g) > K(B) log : 
—r 


for D = 3, say. Together with (9.6) these latter inequalities give 





log »(7, 9) ‘ “ie ® 
ttt. . a | OR a) spa 1 ) 
log M(r,g) a K(p(log~)k o( yj 


> —K(8)loglog M(r,g){1+0(1)}. (9.8) 


Because (9.7) holds for a sequence of values s increasing to 1, (9.8) holds 
for a corresponding sequence of values 7 also increasing to 1, and there- 


fore log p(r q) 
ities: GJ — K(f). 
ae log M(r, g)loglog M(r,qg) . ®) 





Since g(z) is of positive order, the function loglog M(r,g) increases to 
infinity as r—> 1. Hence, if r is sufficiently close to the value 1, 


log u(r, g) < log M(r,g) = —) ae 
log M(r,g)loglog M(r,g) ~ log M(r,g)loglog M(r,g)  loglog M(r,g)’ 





log u(r, g) 
und we have Mest log M(r,g)loglog M(r,g) . 





thus completing the proof of Theorem 1. 


10. A counter-example for Theorem 1 

In order to show that in a sense the result of Theorem 1 is best-possible 
I construct a particular function for which the limit occurring in (1.2) is 
negative. Thus we cannot in general replace the factor loglog M(r,g) by 
any other function which tends to infinity much less rapidly without 
affecting the finiteness of this limit. 
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THEOREM 6. To each finite value B (> 1) there corresponds a function 
m(z) regular and of order B in \z| < 1 such that 





lim sup log u(r, 77) 


10.1 
ns P iog M(r, m)loglog M(r, =) ~ — 


For the proof of Theorem 6, I construct a function regular for |z| < 1 
with zeros of multiplicity [2°"] situated at the points 


a,j (1—2-"+42-29-lyetrik2” (Lo ko ; lon<o). (10.2) 


n,k 


This function is defined as 


on Pp roBn 
x z — l io J 
(2) IIT I] | {Zee ay,)} exp > 2204)" . 
n=1 k=1 m 
m=1 
where p = [8+1] and 
a (a—z) 
a(1—zda) 





Z(z,a) = 1 (10.3) 


[t will be assumed for convenience that 1 < 8 < 2 since the proofs of 
Lemmas 5, 6, 7 need only trivial amendments to apply in the general case. 


Lemma 5. Jf 1 < B < 2, the function 


m(z) = 7,(z) Il II {( l —Z(z,a,, ;.) exp( Z(z, ay) +43Z(z,a,,x)?)}™ 
n=1k=1 








is regular for |\z| <1. 
We must show that this latter product is uniformly convergent in 
each circle |z| < p < 1. Writingz = re anda = |a\ein (10.3) we have 
= 1 +p—ia 
Sic hc ee, (10.4) 
1—za 
When 'z r 1, 'a| < 1, (10.4) gives 
2(1— 
Z(z,a)\ < =( = ) (10.5) 
l—r 
and, when a| > }(3+>p) > #(3+-7), 
we have |Z(z,a)| < 3. We define 


&(z,a) = {1—Z(z, a)}exp{Z(z,a)+4Z(z, a)*}. 
Then, taking the appropriate branch of the logarithm, we have, by a 
Taylor expansion, 
x 1 e 
log &(z,a) < > pa |Z(z, a) mi Z(z,a)\* > Q-m. 
m 


m=1 
m=3 


that is log &(z,a)| < |Z(z,a)), (10.6) 


for jz] =r<p<l, ja} > }(3+ p). 
3695 .2.7 P 
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We can choose a finite value n, such that 


ank Tank > }(3+ p), 
for n > ny and all corresponding values of k. Then it follows from 
0 § 
(10.2), (10.5), (10.6) that, for n, > n, > no, we have 
ne 2" ne 2° ; 9 k 3 
> > [28n) log &(z,a,, ,) < > > 2Bn +3 ny 
n=n, k=1 n=n, k=1 l—p 
2 Ne 
< S > 2(B—2)n 
; (1 —p)s n=", 
For fixed p this latter sum converges uniformly as n, > 00, and therefore 
the product function 7(z) is regular for |z! < p < 1 for each p, which 
proves the lemma. 


11. An upper bound to M(r, z) is stated in 
LemMa 6. The function m(z) satisfies the inequality 
log M(r, 7) < K(B)28%, 


when 1—2-Ytl}=or<l-2-" (N= 


to 
w 
— 


(11.1) 
For fixed @ we write ampa,, , = «,,,, Where |@—«a,,,| < 7, and express 
7(z) as the product of the five factors 
N+1 2" ” 7 
P(z)= TI TT t!—Z4,,)}?" 
n=1 k=1 
P, z) ' II exp[[2°"}{Z(z, a, ,.)+ $Z(z,a,,x)*}] 
. (l<n<N+1; |0—a,,| < 72-), 
P,(z) ~ Il exp[[2°"]{Z(z, a, ,.)+4Z(z,a,4)7}] 
(l<n<N+1; |0—a,,| > 72-*), 
P,(2) — II é(Z,a, ern (n 2 N +2; O—an » < 72 N), 
P.(z) = TJ &(z,a,,)°" (n >N+2; |0—a,,| > 72-%). 


When |z! < 1, \a| < 1, we have, by (10.3), 





| __» |2 2.9 ee “2 1 2 
1— Z(z,a)|? = ‘Attend _ * —#1@ cos(6—«)-+ la ; 
1—2za| 1 —2r|\a\cos(@—«)+r?\a\? 


Hence the modulus of each factor of P,(z) is less than 1, and in particular 


we have log|P,(z)| < 0, (11.2) 
for |z| < 1—2-*. 
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From (10.4) we have 


10 Sts, 0) = pee 
(1—za)(1—Za) 
(1—/\a|){l+rcos(@—a)—r\a cos(6— x)—r? \a)} 
1—2ria cos(@- x)+ ria\? 
2(1 a )(1- ria) 


1 —2rj\a\cos(@—a)+r?\a\2" 





Also re Z(z,a)? : Zz, a) .. a i an — ja )? 
1—2r\a\cos(@—a)+r?\a/? 
which gives 
4(1— \a|)(1—r'a)) 
~ ]—2r\a\cos(@—a)+r7\a 


re{ Z(z,a)+-4Z(z,a)*} (11.3) 





7o° 


In considering this latter expression we shall need later the inequalities 


1 —2r'a!cos(@—a)+r?!a\? (1—r\a})?+-4r|a|sin? 4(@—a) > (1—r{a)})?, 
(11.4) 
and, if |@—a! < 7a, 


| —2r\a\cos(@—a)+r?\a\? > 4rja\sin®? }(@—a) > (27-1)*r|a|(0@—a«)?. 
= (11.5) 
For each factor of P,(z) we have, by (11.3) and (11.4), 
re[28"}{Z(z,a,,.)+$Z(z,a,,)2} < 4.28n, (11.6) 
With each value n (< N+1) there are associated at most two values 
of k which correspond to distinct factors of P,(z). Thus from (11.6) we 


have N41 
log|P,(z)| < > 8.28" < 8, 2AW+2), (11.7) 


n=1 


Next we consider the factors of P,(z) and in particular those for which 
the corresponding value «,, ;, satisfies 


w2-"k’ < |0—a, .| < w2-"(k’+1), (11.8) 
where k’ in an integer and 
l<k’ < 2-1, (11.9) 
By (11.5) and (11.8) we have 
1—2rr, , cos(O—a, ,)+r?r? , > {(0—a,,)/m}? > (k’)?2-2 


since by hypothesis r > }. r,;, > 4. For each value of n and k’ there 
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are at most two values k which satisfy (11.8), and for each such value 
of k we have, by (11.3), 


re[ 28" }{Z(z, a, kb) x $Z(z, a, x)*} < 8. 2B 1 —r,,)(1 =F Tn p)2*™(k’)-* 
< 40, 28"(k’)-2, 


since l—r, . < 2, l—rr,. << 2-*®+2-4t1 < 5,2-", 


Summation with respect to k’ over the range (11.9) and with respect 
to n for 1 < n < N+], leads to 
+1 2%-1 


. N+1 
log|P,(z)i| < ¥ > 80.28(k’)-? < K ¥ 28" < K(p)28, (11.10) 
1k=1 


n n=1 


since > (k’)-? is convergent. 
K=1 
If n > N+2, then, by (10.2) and (11.1), 
Trp > 1—2-4-* > 1—}(1—1r) = }(3+7). 


Hence (10.6) is valid with a = a, ,, and we have 


. (i—r,x\8 one 

log|&(z,a, .)| < |——=} < 2-”. (11.11) 
| i-—r 

To each value n (> N+ 2) there correspond at most 2"-*-!+ 1 distinct 


values a, , for which , , 
: ns 0—«a,,| < 72-5, 
and consequently at most 2”-*-!+-1 distinct factors &(z,a, ,) of P,(z) 
each of multiplicity [2"]. Thus, summing for n > N+2, we have 
x 
log |P,(z)| < p 3 28n(Qu—-N-14 1)2%N—n) 


n= 


to 


x 
92N YY 9(B-2n. 
< 2 - : 
N+2 


R= 


that is log! P,(z)| < K(p)28%. (11.12) 
It remains to examine the function P,(z) and we consider those factors 


for which mk’ 2-N < |0—any| < a(k’+1)2-%, (11.13) 


where k’ is an integer such that 


l<k <2%-1. (11.14) 
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The inequality (10.6) is again valid witha = a, ,,80 that by (10.4), (11.5), 
(11.13), since r > },1r,,, > 4, 
log|é(z,a,,)| < |Z(z,a,%)|* 

8(1 ao rnp) (I a 2rr nk cos(9— Xn) vrs x} 

§ , 2-3"73(9 — hare 

8 . QUN-ayp’)-8; 
that is log |&(z,a,, ;)| << 8.2%>4-™(k’)-3. (11.15) 


For each value n (> N+2) and k’ there are, by (11.13), at most 2"-* 
corresponding values of k, each of which is associated with a factor of 
P.(z) of multiplicity [26"]. Hence by summing over the range (11.14) of 
k’ and for V+2 < n < @ respectively, we have, by (11.15), 
: ¢. ee ‘ 
co 2 a 2 Obn 9AN—n)j]-’\- 
log | P;(z)| < - o 2-2 2 (k’) 


< K 5m, 
since s k-* is bounded. Hence 
log|P,(z)| < K(B)2°%, (11.16) 
and, by (11.2)7(11.7), (11.10), (11.12), we have 
log|z(z)| = log|P, P, P, P, P;| < K(B)28%, 


which verifies the inequality stated in Lemma 6. 


12. As a final preliminary to the proof of Theorem 6 we need an 
upper bound to the function p(r,z). This is provided by 


LemMA 7. The function 7(z) satisfies the inequality 
log p(r, 7) << —2-PN-D log 24-14 K(B)28%, (12.1) 
when r satisfies (11.1). 


We again use the notation of Lemma 6 and indeed several of the 
inequalities which were proved there. However, we must amend 
(11.2). On each circle 


» = Raw 9-N+11L} 9-2N+1 (k =— ] 3. 9N-1__ 1) 
2 2 2 , ; 
there is a zero of x(z) of order [2&%-»]. Hence on each circle |z| = r, 


where r satisfies (11.1), there is some point, say z’, whose distance from 
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at least one of these zeros is not greater than 2-2‘*!, Denoting such 
a zero by a’ we have 


at 


rr. | 
|a-z 


)-2N+1 
enone at woe 
l _ 2’ = 


7\ 


1—2'a’ - 
Since the order of thezero is[2&\-»], the corresponding term of log | P,(z’) 
is less than [2°-’-» log 2-"*+!, and, since all the other terms are negative, 


we have s . ; 
. log | P,(z’)| < 22.N-) log 2N-1+4 log 24-1, 


The inequalities (11.7), (11.10), (11.12), (11.16) are valid with z 2, 
and therefore 


log m(z’) saat DAN Y log 2 ee K(B) 28%, 


which gives (12.1). 
We can now prove Theorem 6 for, by Lemma 6, 

log M(r, z)loglog M(r, 7) < K(B) 2°" log 2%, (12.2) 
for 1—2-*+! < r < 1—2-5, N > 1. Also by (10.2), the number of zeros 
of z(z) in the circle |z| <r < 1 is greater than (1—r)-®-! for r > 3. 
Thus z(z) is of unbounded characteristic in the unit circle since f is 
positive and consequently M(r,7)—- 0 as r—+1. Because of this the 
left-hand side of (12.2) is positive and with (12.1) we have 








___ log u(r, =) _ _ 2A(log 2-1+ K(B)) 
log M(r, z)loglog M(r, 7) ~ K(B) log 2% , 


for 1—2-4*4< r < 1—2-*, N > NN, if N is sufficiently large. Thus as 


N + oc we have that 





: log u(r, 77) , 
lims —K < 0, 
<_ Mr, m)loglog Mr.m) (8) 


the result stated in Theorem 6. 
13. Functions of order unity 


Theorem | is not in general true if we relax the hypothesis that g(z) 
should be of order greater than 1. For example we may consider the 


function ; kt \8 
g,(z) = exp(— =a) (em i) : 
where m is a positive integer, 8 > —2, 
log, x = log(log,,,2) log, x = logz, 


and k is a positive value large enough to ensure the regularity of g,(z) 
for !z} < 1. The function g,(z) is of order unity and it has been shown 
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((3) Theorem 2] that 


l k \B 
log u(r, 91) = — i—r logm 7 — 


_ K(m, B) 7” .. 1 \b 
= 


for r,5 <r <1. Hence it follows that 


log u(r. ¢ 7 l = m l 2 
: I ai By a : — < —K(m, (log ) | | (log, <3) ; 


and the limit which appears in (1.2) has the value —oo. 

The function g,(z) represents an extreme example of the exceptional 
behaviour to be encountered amongst certain functions of order unity 
which possess no zeros in the unit circle. An earlier theorem [(3) 
Theorem 2] can be applied to such functions to give 


- 


THEOREM A. Let g(z) be regular and possess no zeros in the circle \z| <1. 





Then, if g() 1 and 
A m yu 
log M(r,q) - — | | (log, ) I] (og, =) oe <9 < Th 
p=1 u=m 
where each value x,, is real and x,, > 0, there is a constant 
K = K(m, n, fq, %,4)---; X») 
— “ 
such that log u(r,g) > — , 
uch tha og u(r,g) > KA TT ei) 
=m 
jor %,<¢ <i. 


Although Theorem 3 is valid for each function of order unity, the best- 
possible general minimum-modulus inequality which can be obtained 
from it is that for each positive value « and an appropriate constant K, 


log u(r,g) > —K(1—r)-**, 
for a set of values of r which increases to 1. For an improvement of this 
inequality in particular cases a result similar to Theorem A can be ob- 
tained for functions which possess zeros by applying the methods of the 
earlier sections of this paper if we replace (3.1) by the transformation 





. k \8 
l—w = K(1—2)(logm 5 ; , 


where 8 > 0, mis a positive integer and k and K are suitable constants. 
Such amendments do not lead to any particularly interesting general 
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results similar to Theorem 1, and I do not proceed further with this 
aspect of the theory here. 

The constant A(8) which occurs in (1.2) has not been calculated in 
any detail, but a reference to (7.1) and (9.3) shows that the particular 
K(f) is finite for each finite value 8 greater than 1, and tends to infinity 
as 8-1. Although this has not been proved for the function 7(z) 
constructed for the proof of Theorem 6, our brief examination of func- 
tions of order unity suggests that A(8) may not be uniformly bounded 
in general as B > 1. 


I am indebted to Mr. Hayman who pointed out to me the methods 
pursued in the earlier sections of the paper. 
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A NOTE ON THE APPROXIMATE SOLUTION 
OF THE EQUATIONS OF POISSON AND 
LAPLACE BY FINITE DIFFERENCE METHODS 


By J. EVE and H. I. SCOINS (Newcastle upon Tyne) 
[Received 3 April 1956] 


1. Introduction 

One of the numerical methods of finding approximate solutions to 
Poisson’s equation, Laplace’s equation, and to other partial differential 
equations in two dimensions is to replace the continuous independent 
variable by a set of points arranged on a regular grid, usually square, 
and the continuous dependent variable by a function defined on this set 
of grid points. The differential equation is thus replaced by a difference 
equation, which is then solved by relaxation methods. 

If little is known about the solution sought, it is customary to take a 
coarse mesh, solve the finite difference problem for this grid, and then 
decrease the mesh size by a factor 1/Vv2 by introducing grid points at 
the centres of the squares of the old grid. This produces a mesh of 
squares which is orientated at 45° to the axes of the old grid, usually 
parallel to some rectilinear boundaries. It is convenient therefore to 
refer to the old grid as a ‘square’ grid, i.e. parallel to and at right angles 
to the boundaries, and to refer to the new grid as a ‘diagonal’ grid. 

The usual procedure uses the values for the old grid to set up approxi- 
mate values at the new points, so as to give zero residuals at the new 
points { Hartree (4)]; but thereby significant residuals are created at the 
old points. If the finite-difference problem for this ‘diagonal’ grid is 
solved, to a given number of significant figures, by relaxation, it will 
usually be necessary to modify the pivotal values at both the old and 
the new points. The process may now be repeated by decreasing the 
mesh size again by a factor 1/¥2, solving the new difference problem 
on what is now a ‘square’ grid, and so on until no significant change 
takes place in the pivotal values between one stage and the next. 

In following this method numerically for a number of problems, we 
have noticed two significant facts First, for the finer meshes, several 
extra figures have to be carried to allow for rounding errors, which may 
be as large as the number of steps from a given point to the boundary, 


Quart. J. Math. Oxford (2), 7 (1986), 217-23 
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Secondly, if we consider the sequence of values appropriate to a par- 
ticular mesh point at the end of each relaxation, we find that they 
oscillate, those corresponding to a ‘diagonal’ grid forming one sequence, 
while those corresponding to a ‘square’ grid form a second sequence. 

Comment on this type of oscillation in relation to the buckling of 
beams has been made by Salvadori (5), but no explanation is offered. 
We believe we have a partial explanation from a consideration of one 
of the simpler cases. 

We shall show that, if the density function in Poisson’s equation is 
harmonic, the two sequences converge to the solution in the continuum, 
one from above and one from below. Moreover, we suggest, there are 
many cases in which it is unnecessary to follow the usual procedure of 
solving a series of finer and finer meshes until no change is observed in 
the pivotal values. An extrapolation technique can be used to estimate 
the grid size which will be required for a given degree of accuracy, and 
some pivotal values can be accurately determined before relaxation is 
started on a fine grid if the overall solution of the problem is required. 


2. Analysis 
We may write Poisson’s equation in the form 
V2d(x, y)+ p(x, y) = 9, (1) 
where (x,y) is the potential and p(2,y) the density; when p(x, y) is 
identically zero, we have Laplace’s equation. We may set up the corre- 
sponding functional equation 
u(ath,y)+-u(xa—h, y)+ul(x, y+h)+ u(x, y—h)— 

—4u(x, y) +h p(x, y) = 0, (2) 
for a grid with sides parallel to the coordinate axes. Here u(x, y) corre- 
sponds to (x,y) and is a generalization of the point function used in 
the normal difference-equation, being identical with it if (x, y) is restricted 
to the grid points only. Clearly u(x, y) is also dependent on the mesh 
size h. 

For a ‘diagonal’ grid of mesh size h we have the functional equation 
v(ath/ v2, y+h/v2)+v(a+h/ v2, y—h/v2)+0(a—h/ v2, y+h/v2)- 


+-v(a—h/v2, y—h/v2)—4v(x, y) +h? p(x, y) = 0. (3) 


Here v(x, y) corresponds to the potential in equation (1). 

We shall assume that the boundary conditions are that (x, y) is given 
on the boundary, and that this can be applied unambiguously to the 
functional equations. Clearly this will require modification if the grids 
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used do not conform to the boundary, and this can give rise to a diffi- 
culty which we shall discuss later. 

Little seems to be known from the pure-mathematical standpoint 
about the nature of such functional equations in two dimensions, so we 
assume that solutions exist which are regular functions of x, y, h for 
values of / sufficiently small. Under these cireumstances we can expand, 
where necessary, using a Taylor series and write, following the notation 


of Woods (8) and Bickley (1). 


h* h4 ‘ 
| v? ra 29") sagt 374V?)4 O(n) | u(x. ele.) - 0, (4) 
- odD 
and 
| v2 hv 174) 4 an (V8+1294V?) + O(A®) |u(a, y) + pla, y) 0 
24 1440 ii . 
(5) 
where V2 Ae n ae G* a (6) 
Ox* —soy* CXCEY 


The functional difference-equations (2) and (3) have thus been replaced 
by differential equations, similar to the original differential equation 
(1), but involving correction terms which differ in the two cases of 
the ‘square’ and ‘diagonal’ grids. As has been mentioned by one of us 
elsewhere (2), there is a danger that solutions of differential equations 
of degree higher than the second may creep into the solution here, but 
this has yet to be observed in practice for what is essentially a jury 
problem. 

Since only h? is really involved in the equations, we may suppose that 
u(a,y) and v(2,y) can be expanded in a power series in h*, say 

u(r, Yy) = do(X, y) +h? db, (x, y)+htG.(%, y) +... | (7) 

v(x, Yy) W(x, y) + hh, (x,y) h+p.(x, y)-+ ~ \” ; 

Substituting into the differential equations (4) and (5) and comparing 

coefficients of powers of h?, we see immediately that ¢o(z, y) and &o(x, y) 

must both be the required solution ¢(x, y) of the original equation (1). 
Moreover, ¢,(2,y) and y,(2,y) are given by 

V2d, (2, y) = —3h(VI—274)d(a, y) | 

V2b, (x,y) = —3,(V1+474) (2, y) \ 

respectively, where the boundary condition on both ¢, and ¥, is that 

they vanish on the boundary, leaving no ambiguity as to the solutions 

required. In certain special circumstances these principal correction 

terms are equal and opposite; they are, at the moment, corrections to 


(8) 
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the solution of (1) such that equations (2) and (3) are satisfied. This is 
the case when V*p(x,y) = 0, implying V4¢d = 0. This is, of course, also 
true for Laplace’s equation, i.e. when p = 0. The second-correction 
terms (i.e. the h* terms) can be shown to be approximately equal when 
the density is harmonic, and equal when the density is identically zero. 

From this we are led to suppose that, when the finite-difference 
equations corresponding to the functional equations (2) and (3) are 
solved with the correct insertion of the boundary conditions, then the 
pivotal values w; at a particular point (denoted by 7), considered as a 
function of the mesh size h I /n, can be approximately represented by 


w (rn) A, +B, h?+0C,h'+..., (9) 


where the coefticients depend on the particular lattice site and the alter- 
native signs are taken positive for ‘square’ grids and negative for 
‘diagonal’ grids. 

In the case of Laplace’s equation this hypothesis is correct, and, when 
the density function is harmonic, at least the first two terms are correct 
for Poisson’s equation. For the general density function, the usefulness 
of the hypothesis will depend on the relative magnitudes of Vp and 
744, as can be seen from equation (8). 

The order of magnitude of B; can be estimated by an approximate 
integration of equation (8) using data obtained from a coarse mesh, and 
this can be used, if required, to estimate the value of the mesh size / 
necessary for the difference between the finite-difference solution and 
the continuum solution to be negligible. 

The validity of the hypothesis embodied in (9) is subject to the fairly 
severe restriction that the boundary has in some way to conform to the 
grid. There must be no ambiguity as to how the boundary conditions 
imposed on the potential d(x, y) in (1) are to be applied to the difference 
equations, as distinct from the functional equations (2), (3). It is 
preferable that the boundary points of the grid should lie on the 
boundary curve. 

That the restriction is very necessary is shown by considering the 
special case where the boundary is a circle. Clearly, in this case, there 


is no distinction between ‘square’ and ‘diagonal’ grids, so that no 
oscillation should be observed. However, the way in which the boundary 
conditions are applied to the finite-difference equations will modify the 
problem, and variations in the solutions as h? is decreased will depend 
on the way the boundary conditions are applied and not on changes due 
to grid orientation. 
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Thus we must restrict attention, for the strict application of (9), to 

cases where the oscillations of the solutions, as h? decreases, are due 

solely to changes in the grid orientation, and not to changes in the 
manner in which the boundary conditions are applied. 


3. Example 

The problem of the solution of Poisson’s equation inside a square for 
constant density and zero boundary value has been considered by several 
authors, (3), (8), and an algebraic expression for the potential is known, 
(6), (7). 

We have considered the problem of finding the potential at the centre 
of the square to an accuracy of 1 part in 10* following the methods 
discussed above. 

Taking the differential equation |cf. Hartree (4) Ex. 37] in the form 

V2V 2/a?, 
inside a square of side a, and V = 0 on the boundary of the square, the 
difference equation takes the form 


“wen ’ » 2 
= we; 4u,-+-2n 0. 
) 


Here w; represents a typical pivotal value corresponding to either type 
of grid, S w, denotes the sum of values at the four immediate neighbours 


of i, and n = ah~', h being the mesh size. 

Since we wish, in this instance, merely to examine the values of w 
corresponding to the centre of the square, only a restricted class of grids 
need be considered. For ‘square’ grids, n* is the number of small squares 
into which the basic square is divided, and m is given only the values 
2, 4, 6, 8,.... For ‘diagonal’ grids, n/v2 1, 2, 3, 4,... give grid points 
at the centre. 

The values of w corresponding to the centre of the square, and the 


various values of n, are shown in Table 1. 


TABLE 1 


2 4 6 8 
0-12500 0-14062 0-14425 0-14560 
s2 2v2 3.2 4\2 
w 0-25 0- 16667 0-15476 0-15127 
For n 4,2, 6, 8 the difference equations were solved by relaxation, 


the other cases being solved analytically. 
The extrapolation technique used elsewhere by one of us (2) for 
extrapolating the eigenvalues of sets of difference equations to the case 
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of zero mesh size can be used here for an extrapolation to h? = 0 based 
on equation (9). 
If w(m) and w(n) correspond to n-values m and n, then the ‘linear’ 
extrapolant w(m,n) is given by 
wim, n) = sa ee 
n°-—- m= 
the negative sign being taken if m and x correspond to the same grid 
type, and the positive sign being taken if m and n correspond to different 
tvpes of grid. 

This modification of a Neville extrapolation |(4) 85] can be extended 
quite straightforwardly to quadratic and cubie extrapolation, and a 
Neville table can be set up to display the extrapolants. Table 2 shows 
the extrapolants for the natural sequence of grids described in our 
introduction, while Table 3 shows the extrapolation for the four finest 


erids considered. 


TABLE 2 


) wn) wim, n) wl, m,n) 
v2 0-25000 
-16667 
2 0-12500 “14815 
-15278 
3y2 0- 16667 “14814 
-14930 
4 0-14062 “14719 
-14772 
4v2 0-15127 “14741 
-14749 
Ss 0-14560 


TABLE 3 


n w(n) w(m, n) w(l,m,n) w(k,l, m,n) 
3v2 0-15476 
*14775 
6 0-14425 14729 
*14755 -14738 
4,2 0-15127 -14741 
-14749 
8 0-14560 


The relaxation process leaves a significant rounding error, which in 
the case of n = 8 may be as much as 4 in the last significant figure, and, 
since the extrapolation technique tends to emphasize such small errors, 
we draw the conclusion from these two Neville tables that, in the 


continuum, the potential at the centre is 0-1474+.0-0001, which is in 
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agreement with the value 0-14734 obtained from the analytic solution 
and the value 0-1473 found by Fox (3) (after sealing). 

The extrapolation technique applies equally well to other grid points 
which appear in several grids, the ‘square’ grids giving lower and 


‘diagonal’ grids higher values than that corresponding to the continuum. 


4. Conclusion 

We have shown that the pivotal values at a particular point may be 

approximately represented by 

w = A+ Bh?+-O(h4), 
where A and 6 depend on the position of the point and the alternative 
sign is taken when a ‘diagonal’ mesh is used, as compared with a ‘square’ 
mesh. 

The approximation will be perfectly satisfactory if the density function 
p is harmonic, and will still be useful if V*p is small compared with 74¢. 
The main restrictive feature is that the change from ‘square’ to ‘diagonal’ 
grid should not affect the way in which the boundary conditions are 
applied, and so will have straightforward application only for problems 
with rectangular boundaries. 

The advantage of the method is that in a particular problem the value 
of the potential at a few points can be fixed with some accuracy, and 
these can be used to obtain a more detailed picture of the solution using 
either a fine mesh, or the more sophisticated methods of Fox (3). In 
this way many of the difficulties of the relaxation process for fine meshes 
can be avoided, especially if computing facilities are available for solving 
accurately the finite-difference equations for the smaller values of h, 
where the number of grid points is not more than (say) 40. 
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FOURIER SERIES WITH GAPS 


By P. B. KENNEDY (Cork) 
[Received 9 May 1956] 


1. Let {n,} (k = 1, 2, 3,...) be a strictly increasing sequence of positive 
integers. In this note f(x) will denote a real function, L-integrable over 
(—7,7) and having period 27, whose Fourier coefficients a,,, b,, vanish 
except possibly when x has any of the values n,. Noble (1) supposed 
that gatas oe 
Meri Mk on (k-> 00), (1.1) 
log n;. 
and derived a number of properties of the coefficients a,, 6, under 
various hypotheses about the behaviour of f(x) in an arbitrarily small 
interval. Noble’s proofs start from the fact that 


7 


a7 


a ie : 
a, =— | f(x) P(x) cosn,x dx 


— 


(with a similar formula for 6,,,), for any trigonometrical polynomial P(x), 
with constant term 1, of degree less than min(n,—7,_,,%j..,;—”,). He 
shows that, given «a and 8 with —a7 < a < B < z, one can choose P(x) 
so as to be very small in (—z, x) and (8,7); |a,,| is thus made to depend 
mainly on the behaviour of f(x) P(x) in (a, f). 

More powerful methods of approach to this kind of problem were 
developed by Paley and Wiener (2). In this note I show how these 
methods can be used to give quite simple proofs of Noble’s results, 
under the less restrictive gap hypothesis 


Nii —Ny, > 00 (k->00). (1.2) 


2. Following Paley and Wiener (2), we consider, in the first instance, 
a more general situation than that described in § 1. Let 
fA,} (—20 <k< a) 


be a sequence of real numbers satisfying 


Anay >A, (—0O<k < ~w), Apat—Ayz > 0 (\k} +00), (2.1) 


Quart. J. Math. Oxford (2), 7 (1956), 224-30. 
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and let {A,} (—oo < k < ) be a sequence of complex numbers such 
that . 

> |A, |r <0 (O<r< 1). (2.2) 
Put d(r,x) = > A, reir (0 <r <1) (2.3) 


for all real x; the existence of 4(r,2) is ensured by (: 


to 
to 
— 


Suppose now that 


d(x) = L*-lim¢(r,x) (\w—ay! < 8), (2.4) 
r—1 
where 2 is fixed and 6 > 0. We then have the following theorems. 
THEOREM I. Suppose that d(x) has bounded variation in \~—2,| < 5. 
Then 


A, = O(\A,|-2) (\k| > 00). 


THEOREM II. Suppose that d(x) € Lipa in |ja—ax, 5, where 


0 
O << ad. 
Then A, = O(\A,|-*)  ({k| > 00). 
THEOREM III. Suppose that d(x) € Lipa in |~—a,_| < 8, where 
t<a<l. 
Then P A, < ©. (2.5) 
THEOREMIYV. Suppose that, in \x—ax, < 8, d(x) has bounded variation 
and d(x) € Lipa, where 0 <a <1. Then (2.5) is true. 


Returning to the situation described in § 1, we deduce from Theorems 
I-IV the following theorem on Fourier series with gaps, which comprises 
Theorems 1, 2, 3, and 5 of (1), with (1.1) replaced by (1.2). 

THEOREM V. Let {n,} and f(x) be as in§ 1, and suppose that (1.2) holds. 
Then the following statements are true. 

(i) If f(x) has bounded variation in some interval I, then 
a,, 0 O(n-") (n>). 

(ii) Zf f(x) € Lipa in some interval I, where 0 < «a < 1, then 
a,, 6, = O(n) (n> 0). 


(iii) If f(x) € Lip« in some interval I, where } <<a <1, then 


(\a,,|+|b,|) < oo. (2.6) 


n 


=hA8 


(iv) If, in some interval I, f(x) has bounded variation and f(x) € Lipa, 
where 0 < a < 1, then (2.6) is true. 
3695.2.7 Q 
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Without the hypothesis (1.2), and with (—z,7z) for Z, Theorem V is 
classical [(3), 18, 135]. 

3. The key to the proofs of Theorems I-IV is the following lemma, 
which is essentially due to Paley and Wiener (2). 


Lemma 1. Let {A,\; {A,}, (x), a, and 8 be as in § 2, but let (2.1) be 


re place d by A, +47 r;,. > 87d 1 ( -cO <. k . 76 : (3. l ) 
Then SA,2< 88 [  [d(x) Pde. 
: r—2o!<8 


To prove Lemma I, put 
> K 
d(r,a;K) = > Apr eer (0 <r < 1). (3.2) 
“k 


Then (3.1) and the argument by which Paley and Wiener proved the 
inequality (31.10) of (2) lead to 

g -, 9 sin*{475-1(2—2)} a ; rae 
d(r, a; K)|2— - rods dx > 278-1 > |A,|?r?e, 
ae (7—2y)? Kk 


r—Xl <0 
and this, with the inequality 
a d*sin*{ 4775 1(x—2y)} 


l PETS AL. 


~ 167?(x—2)? 


, 


gives ( d(r, x; K)|\* dx > 36 s A, |?r2s, (3.3) 
r—%\<8 iia 

But it follows from (2.2), (2.3), (3.2) that, for any fixed r such that 

<7 4% d(r,x;K)>d(r,x) (K >) 

uniformly in every finite range of values of x. We therefore deduce 

from (3.3) that 


x 
| d(r,x)|2 dx > 385 > |A,|2r2%, 
xr—ao| <8 i 


Lemma | now follows since (2.4) implies 


[  jd(r,x)itdxe> | |d(x)@dx (r>1). 
r—%9|\<8 r—%9\<8 

The deduction of Theorems I-IV from Lemma 1 proceeds by fairly 

straightforward generalization of the arguments used to prove the 

classical analogue of Theorem V; the part played by Bessel’s inequality 

in the classical theory is played in our proofs by Lemma 1. We assume 
throughout that 

Ania —Ay > 1675-1 (—00 < k < cw). (3.4) 














FOURIER SERIES WITH GAPS 227 
In view of (2.1), this can be achieved, if necessary, by adding to (zx) 
a polynomial in e«", a process which affects neither the hypotheses 
nor the conclusions of the theorems. 


From Lemma | we deduce the following lemma, which we need for 


the proof of Theorems IT, III, and IV. 


LeMMA 2. Let {A,}, {A;,}, d(@), 2%, and 8 be as in § 2, but let (2.1) be 
h k 0 3 
replaced by (3.4). Let the integers j, K satisfy 


Ax | > 276-1, 0 <j < 4} 3/Ax). (3.5) 
2) | 2)—1)z 
Put Ox) bla J db) ae (5) \, (3.6) 
Ax |) t Ax! J 
Then S,. > A,|?< 88"? | p(x) * dx. (3.7) 
b Agi <iAkls Ag x— Io! <4 
To prove Lemma 2, put 
: 1(2)—4)zA,) . aA, 
A,()) 21 A, exp!" J~3) KY sin _ kK (-a<k<o). (3.8) 
\ AK 2 \AK 
Since |A,(j)| < 2.A;!, we have, by (2.2), 
Y |A,(j)\r* <0 (0<r<]l) (3.9) 
Put (7, x) > A,(jrrerX (0O<r <1). (3.10) 


The identity 

wah, blr, a (j—1)n) 

Axi) "| Aci J 

is easily verified, and from it, combined with (2.4) and (3.6), we obtain 


bx) L*-lim¢,(r,2) (|w—a9| < 48), (3.11) 


ok 


since, by (3.5), the points 


x+(2j)—1)z|A, |=}, x+2)7 AK 1 
ie i interval (2 —5,x%+5) whenever |x—2,| < 46. It ’ follows 
lie in the interval (2)—8,.x7)+5) whenever |\x—a,| < $8. now follows 


from (3.4), (3.9), (3.10), (3.11) and Lemma 1 (with A,, 4, 5 replaced 
by A,()), d;, $6) that 

> |A;,(j))? < 168-1 f;(x) ? dx. (3.12) 
' ? IIo <5 


But by (3.8) and the obvious inequality 


7A k ~ 1 ( 1 


— 
sin* 
» 


- AK 
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we have > |A,(9)|? > 28x, 


where S, is as in (3.7). This, with (3.12), gives (3.7), and so Lemma 2 


is proved. 


4. Proof of Theorems I-IV 
Throughout this section the integer A is assumed to satisfy the first 
inequality of (3.5), and the symbol ¥ denotes summation over all 
j 
integers ) satisfying the second inequality of (3.5). 
For the proof of Theorem I we do not use Lemma 2, but employ the 


notations introduced in that lemma and in its proof. Put 


At = > A,(j) (—0o <k < ow). (4.1) 
j 
Then y Atirli< ao (0<r< 1). (4.2) 
Put $*(r,x2) = > Af rrvede = ¥ 4,(r, 2), (4.3) 
=x 7 
so that, by (3.11), 
L?-lim d*(r, x) = ¥ $,(x) = $*(x), say, (4.4) 
r—1 j 
over the range |x—2| < 48. It follows from (3.4), (4.2), (4.3), (4.4), and 
Lemma 1 (with A,, 4, 5 replaced by Af, ¢*, $5) that 
At)? <3 At? < 168-2 [|p *(x) |? de. (4.5) 
: xz Tol 45 


From (3.8) we have 


A,(j) = 24g (0 <j < }18|Ag)), 


and so, by (4.1), A |? > 4-78? Ag |? | Ag |?. (4.6) 
Also, for x—x,| < 48 and integers j, A satisfying (3.5), the intervals 


Ll yi} 95 


fx (2j)—l)a Ax 


are non-overlapping subintervals of (2)—8, 2 )+5). Therefore, if V is the 


total variation of 4(x) in |x—2,| < 5, we have 
p*(x)| < ¥ |b(x)| < Vi (|w—aXq|< 48). 
J 


This, combined with (4.5) and (4.6), gives 


2 — O(\Ax|-2)  (|K| 00). 


Ax 


This proves Theorem I. 
To prove Theorem II, we apply Lemma 2 with any fixed integer j 
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satisfying (3.5). We observe that S, > |A,-|*?, while the right-hand side 


where w(t) is the modulus of continuity of 4() in (#)—8, 2 )+8). There- 


of (3.7) is at most equal to 


fore, by the Lipschitz condition, 
Ay |? = O(Ag|-*%) (| K| > 0). 
This proves Theorem II. 
To prove Theorem III we start with 
S,. O(\Ax|-**) (|K| > 0), 
which we obtain from (3.7) just as in the preceding paragraph. Let p 
be a positive integer. Either the set of integers & for which 


2? < |A,| < 27+! 
is empty, or there is a member of this set, say K = K(p), which has 
largest modulus; and in the latter case the set is included in the set of k 
for which }/A,;!| < |A,! < |Ax!. Thus, in either case, 
oz Ay? < Sx = O(/AK|-™*) = O(2-*”). 
— 


From this we obtain (2.5) by a familiar application of Cauchy’s inequality. 
In fact, by (3.4), the number of terms in the summation is O(2”), and so, 


by Cauchy’s inequality, 


. | pp Vv 2 }) fou x)p 
> 4g = Ol” Y 14, )!) = Of20-a}, 
2 Mt) <2 2P<\Ag)< 277? 
Now (2.5) follows at once, since a > $. This proves Theorem IIT. 


To prove Theorem IV we again apply Lemma 2. By the hypotheses 
of Theorem IV and the sentence immediately following (4.6), we have 


> |6;(x)? s Veo O(\Ax|-*), 
AK 


\ 


uniformly for x—.2,' < 38. Therefore, summing (3.7) over 
YaTa 
0 <j < jr Ag 
we obtain hor 1§ AK Sk O( AK a 
and so S,; = O(|\A;,\~!-*). We can now deduce (2.5) from this by arguing 


with Cauchy’s inequality, as in the proof of Theorem III. This proves 


Theorem IV. 
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5. Proof of Theorem V 
To prove Theorem V, let {n,}, f(x), a,, and b, be as in§ 1. Put 
A, 0, A, (ay, ib,,) (i 4 0), A, A k (k ~ 0), 


and A. =n, (k > 0), A. A, (k< 0). 


Then (2.1) and (2.2) are satisfied, and 
d(r,x) = ¥ (a, cosnxr+b, sinner". 
i 


Theorem V now follows at once from Theorems I-LV, with 4(2) = f(x), 
by virtue of the following simple lemma. 
Lemma 3. Jf f(x) € L°(1) and x—2x,y| < d is a proper subinterval of I, 


the . 
‘ | < 8). 


f(x) L?-limd(r,x)) (a—x 


r—1 

This lemma is almost certainly known, but I indicate the proof since 

I cannot quote a source. 

We suppose, as we may without loss of generality, that J is a sub- 
. 
interval of (—7z,7). Put 
fi(z) = f(x) (xe TL), f(x) 0 otherwise, 

Jo(x) = f(x) —fi (2). 


Let a, BD and a, b! be the Fourier coefficients of f,(7) and f,(x) 


vt it 
respectively, and, for j = 1, 2, let 
AM (7, 2) tay’ 4 > (a? cosnx+b” sinna)r® (0 <r <1). 
n=1 


Then, since f(a) ¢ L?(—7z, 7), it follows from a known theorem [(3) 87] 

that f(z) = L?-lim g(r, xz) (|x| < 2). (5.1) 
rl 

Further, by Fejér’s theorem on A-summability of a Fourier series 

| 5 re have 2 . . 

[(3) 51], we have pr, x) > f(z) (r>D), (5.2) 
5. Since 


A(r,2) = dr, x)+4°(r, x), 


Lemma 3 follows from (5.1) and (5.2). This proves Theorem V. 


uniformly for «—2x, 


I wish to thank the referee for his suggestions. 
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ON AN EXPANSION IN EXPONENTIAL 
SERIES 


By 8. VERBLUNSKY (Belfast) 


[Received 16 May 1956] 


1. Introduction 
Let k(w) be a function of bounded variation, possibly complex, defined 
in the interval (0,1). We shall be concerned with the expansion of a 
function in a series of exponentials e, where the A, are the zeros of the 
integral function 1 
A(z) = | k(w)e™ du. (1) 
0 
We may suppose that 
k(0) = k(0-+), k(1) = k(1—). (2) 


The method of investigation is based on the theory of residues, and 
requires some information about the order of magnitude of A(z) on the 
path of integration. This is attained by restricting k(w) to satisfy the 


condition 


k(1)k(0) 4 0. (3) 
[ shall suppose that all the zeros of A(z) are simple, i.e. that 

A’'(A,) = 0, (4) 
and also that A(0) = 0. (5) 


Conditions (4) and (5) could be relaxed without difficulty. The formulae 
would then become complicated. Choose the notation so that 


0 “.. A, . A, 
k ' : 
Let d,(u) = ku), ,(u) = - (w) Au | k(v)e” di (6) 
%,(u) = 1, p(w) = ed, (7) 
wherev = 1, 2,.... Then the functions 4, ys are bi-orthogonal in (0, 1), and 
: ; 
| doo du = A(0), db, ws, du = A’(A,). (8) 
0 0 


Quart. J. Math. Oxford (2), 7 (1956), 231-40. 
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For a function f(u) ¢ L(0,1) we write 


1 
] - 
i= AO | f(wWdo(u) du, (9) 
0 
com 
y= FL | f(w)d,(u) du (v = 1, 2,...). (10) 
0 

Let B(z) = zA(z). (11) 

: 
Then B(z) = k(1)e*=—k(0) —*| e*™ dk(u). (12) 

0 


By (2) and (3) it follows that there is a positive constant C such that 


1 | h(1)e? p> ©}, 
B(z)| > , sm (13) 
L4)k(0)|  (w@ <—C), 
where z = x+iy. The zeros of B(z) are therefore in the strip x) < C, 


and a simple argument [see Lemma 2, § 2] shows that the number of 


zeros of B(z) in the strip 
2n7 <y < (2n+-2)z 


is bounded with respect to the integer x. We can therefore find a positive 
number 8 such that, if each zero of B(z) is the centre of a disk of radius 5, 
then there is an unbounded increasing sequence of positive numbers r,, 
such that the circle C,,: |z| = r,, has no points in common with the disks. 
There is a wide choice for the r,: e.g. they can be chosen so that 
n.14—lyp <1. Let v, denote the greatest integer v such that |A,| < r,. 
I shall prove 


THEOREM 1. Let f(u) € L(0,1). Then, as p>, 


vp 3 sin r,,(t—w) 


to + Ya,e*—— | f(u)—* du 
1 





7 —1 


converges to zero uniformly in any closed interval (x, 8) interior to the open 
interval (0,1). 
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An easy consequence is 


THEOREM 2. Let f(u) be integrable in every finite interval, and satisfy 


the equation , 
| k(u)f(t+-u) du = 0 (14) 
0 
for all t. Let (a,b) be an assigned finite interval. Then, as p +>, 
b 
; 1 f ,, ,sinr,(t— 
> a, e*— | f(u) ~~ du 
I me t—u 
a 


converges to zero uniformly in any closed interval interior to the open 
interval (a,b). 

A slight change in the proof (see § 5) shows that, if the set of distances 
between pairs of zeros of A(z) has a positive lower bound, then we can 
put v, = p in Theorems | and 2, and replace r, by |A,]|. 

The case of Theorem | in which f(w) is of bounded variation in (0, 1) 
and satisfies the condition 

1 

| A(u)f(u) du 0, 

0 
as also the case of Theorem 2 in which f(x) is of bounded variation in 
every finite interval, has been considered by Delsarte (1). It should be 
remarked that Delsarte subjects (uw) to a further restriction, viz. that 
the continuous part of (wv) shall be absolutely continuous. The proof 
of Delsarte is rather long, and would seem to be incomplete. [See (1) 438, 


line 15.] 


2. The function B(z) 

It will be convenient to introduce the following 

DEFINITION. An integral function H(z) will be said to have the ‘d-n 
property’ ona set S if, given 8 > 0, there is a positive number y (= 7(8)) 
such that, if each zero of H(z) is the centre of a disk of radius 5, the n, for 
all ra in NS outside the disks. H(z) P n- 

If S is the sum of sets S; (j Rocks n) on each of which H(z) has the 
5-n property (with different functions »(3)), then H(z) has the 8-y 
property on NS. 

LemMA 1. B(z) has the 8-n property on the whole plane. 

Proof. By (13), B(z) has the 5-y property on the set x > C and on 
the set 2 : C. It suffices to prove that B(z) has the 5-7 property 
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on the set |x| < C. If not, there is a 5 > 0 and a sequence of points 


C, (= €,+77,), with |€,| < C such that 
‘ Tha hm - 
1, > 0, B(f,) > V, B(z) =U 
for |z—{,| < dand |x C. The condition |x| < C may here be omitted 
since b(z) = 0 for |x| > C. We may suppose that 7, -> 00 or —00; say 


the former. 


Let #,, denote the rectangle 
2) = 0, 2mm < y¥ < (2m+-2)r (m = 9, 1.,...). 
To each v there corresponds an m, such that ¢, ¢ R,,,. Write 
a C,—2mim,. 
We may suppose that the sequence { is convergent; say, ¢, > ¢’. Let 
tha l 
F(z) = B(z+-2im) = k(1)e*=—h(0) — | eee? dk(u). 
0 
Write G,(z) = F,,(z). Then 
G,(¢,) > 0, (15) 
Gz) ~0 for |z—{| <6. (16) 


Consider the functions G,(z). They are uniformly bounded in the 


rectangle ’ 
R*: |x| < C+, —l <y < 27+1. 


Let RF’ denote the rectangle 
r| < C+4, —k<cy < 2n+4H. 


There is a sub-sequence G, (z) and a function G(z) such that 


(i) G,(z) > G(z) uniformly in R’; 
(ii) GY?(z) > Gz) (n = 1, 2,...) uniformly in R’, 


By (15) and (i), G(¢’) = 0. Now G(z) = 0. For 
; 
G™)(z) = |k(1)e? = yrerte2rimu dk(u)| (m = m,,) 
‘ 0 ' 


> 4/k(1) |e-© 
for |x| < C if n is sufficiently large. It follows that there is a zero of 
G, (z) in any assigned neighbourhood of 2’ if r is sufficiently large. This 
contradicts (16). 
Lema 2. The number of zeros of B(z) in R,, (n = 0, +1.,,...) ts bounded 
with respect to n. 


Proof. If not, there is a sequence m, (v = 1, 2,...) with |m,| > oo such 
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that the number of zeros of B(z) in R,,, tends to 00. We may suppose 
that m,, — 00 or —oo, say the former. Write 

G, = F, T= & 


my 


typ 


os 
There is a sub-sequence G,, tending uniformly in R’ to G (+ 0). Hence 
the number of zeros of G,, in Ry (i.e. the number of zeros of B in T,) 
does not exceed the number of zeros of G in R’ if r is sufficiently large. 
This is a contradiction. 

The contour C,, of § 1 will be traversed in the positive direction. The 
parts for which x > 0, x < 0 will be denoted by C3, Cy respectively. 
By Lemma | and (13), we have, as p > 0, 


Be = OW on Css (17) 
7 si 
Bo O(1) on Cz. (18) 
3. Proof of Theorem 1 
1 
Write g(t) | k(u)f(t+-u) du. (19) 
0 
Consider the expression 
1 u 
1 ett : -; 
I ee | A dz | (ue du | f(vje* dv 
Cp 0 0 
— 2 ; 
se | am % | f(vje** dv | k(uje™ du 
Cp 0 v 
. Aut : ] )) 
— en - ; L(y ' 
— db, (v)- di 
z. 700 | Heyer) — 
J 0 
~ Apt (0) XS ed! (20) 
e , . —~ — 
qe NLA) 
We can write ; 
ntl, . T 
C; ct 
= Ji +d). (21) 
. 2 
Now z| k(uje™ du = k(1)e*—k(v)e” — | e™ dk(u). 


v v 
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Hence 
; 
=| BO dz je fl =e —k(v)e? tf ean) dv 
Cr v 
(1) 
- k(1) e i: in a dy — 0 | +. —_— 
=| Bo): ds | F (eye dv —g(0) Bo“ 
¢, c; 
1 u 
;_€ dz [ dk(u) [ fivjexe—” dv 
oa > etu-v) dy 
| ae | roy 
C5 0 0 
= K,—K,— ky. 
We recall that an integral of the form 
q 
| d(u+p)e™ du, 
p 


where yu, p, g are bounded parameters such that gq > p > 0 and ¢ is 
L-integrable in the interval (min p+, max q+), tends to zero uniformly 


as |z! > oo with x < 0. Further, 


| |e* dz| = O(1) 
C5 
as p—>oo if A is a positive constant, and is uniformly O(1) if A is a 
parameter bounded below by a positive constant. 
Consider A,. On writing 1—v = w, we see that 
; 
| f(v)e 1-0) dy 
0 
is uniformly 0(1) on C>. Nowt >a > 0. By (18), A, is uniformly 0 (1). 
Consider A;. Then 
u 
M(u) = | f(eje-? dv 
0 
is uniformly 0(1) on C>. So then is 
1 
| M(u) dk(u). 


0 


Hence, as before, A; is uniformly o(1). 
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. lag A(z) [10 


Cc 


i v 


0 0 


1 


[ et dz [ fiw 0 dy — [ i 


. . - 
( 0 C; 


howd. 


Then 
1 


L, f(v) dv | e*t-v) dz = 
0 C, 
. [ k(u)e™ du 


k(v)e= 


Since 


we have 
1 


los | B(z) 
C; 0 

- oct 

™ 
q(0) —— az 
B 


where 
~ | KO 55 a 
| — ; 
7) Biz) ; 
( 0 
In 1, | fe 
0 


is uniformly o(1)on C}. Now 1—t > 


o(1). In JL, 


u 


is uniformly o(1) on C. So then is 


1 
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27 [ f(r) 


dz | f (ve “f(o}e*—K(0) — 


= M, —M,, 


= 1—Bp a 0 


1 
)= | flee 
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| } _ [ |x “ dudv 


v 


1 
dz | floje-* | k(u)e™ dudv 
. 4 


0 


1 

sin r,(t—v) dv. 
t—v 

0 


v_ (0) — ( e" dk(u), 


0 


v 


, ezu aw) dv 


. 


0 


z | f(v)e-*" dv, 


. By (17), Mis uniformly 


2(u—v) dv 


( N(u) dk(u). 


0 
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Hence .V, is uniformly o(1). Thus 


ext j 
j BR te,(t), 


Cp 


J, +d, = L,—g(0) 


where ¢,,(¢) is uniformly 0(1). Since 


l [ zt le ] XS eArvt 
: . a 2 ie 
2m J Biz) A(o) ' ~ A, A’(A,) 


Cp 


and a») = g(0)/A(0), this proves the theorem. 
With the notation 
: : 
f(we-** du | k(v)e*" dv, 


. 


u 


* 


l 
y= 
. Ae 
0 


we have the following 


CoroLuaRy. Under the conditions of Theorem 1, if g(0) = 0, then 


1 
Vp lf, .sinr,(t—u 
Fate} f py) Bol—9 gy 
: eae t—u 


. 


0 
converges to zero uniformly in (x, B). 


For a) = 0 and g(0) 
x —a* = — - = 0. 


woo FS 
4. Proof of Theorem 2 
We have g(t) = 0 for all t. In particular, g(0) = 0, so that a, = af. 
It suffices to prove that, if y is any assigned number, then 





y+ 
>. 1 ff sin 7, (t—u 
> at et _ — flu) a rolt—) du 
i w t—u 
Y 
* 


converges to zero uniformly in any closed interval (y+, y+) interior 
to the open interval (y,y+1). Write uw = v+y, f(u) = F(v), and 


1 1 
pt = 7 | F(v)e-*? dv k(uje?" du. 
0 v 
Since | k(w)F(w) dw = g(y) = 0, 


0 
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we can apply the corollary to the function F(t). It follows that 


1 
Up Lf —,..smr -v 
¥ pe dt | F(v) 0 *pt— 2) dv 
I 7 atlas 
0 


converges to zero uniformly in («, 8). Hence 


2 lf. .sinr,(t -u) 
» B* Avy eAvt _ f(w) oh eens SF 
| u 


converges to zero uniformly in (a+y,8+y). It remains to prove that 
yk \ 7 
VD, é ty 


For a fixed v, consider the absolutely continuous function of w, 


1 w+u 
D(w) | k(u)er" du | f(s) Avs ds, 
0 w 
1 
Then D(y) | he(w)er du | Fi(v)e AMe+y) dy 
0 0 
A’(A,) Bee Avy, 
and @(0) = A’(A,)a*. 
On writing 
1 h 
l F re 
[O(w+h)—O(w)] k(wye ~s | f(wtut+s)es dsdu — 
f} 7 
0 0 


1 w+th 
i if ce 
we ki u ye Apu du i f(s)e-*’s as, 
(} 


. 


5 0 w 
we see that, for almost all w, 


Y'(w) ert g(w)—f(w)A(A,)] = 0. 


Thus O(y) (0). This proves Theorem 2. 


5. Suppose that the set of distances between pairs of zeros of A(z) 
has a positive lower bound d. Let 8 < 3d. Let each zero of A(z) be the 
centre of a disk of radius 8. We denote the disk whose centre is A, by 
A,, and its perimeter by y,. We define the contour C’, to be the circle 
fi A,,!, modified as follows. If I, meets A,, where v < n, we 
replace the are of ., which is inside y, by the are of y, which is outside I. 
If ., meets A, where v > n, we replace the are of I, which is inside y, 
by the are of y, which is inside [,. Every point of C), is at a distance 


a 


not less than 5 from all the zeros of B(z). Hence, as n > 00, e*/B(z) 1s 
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O(1) on the right half of C’,; and 1/B(z) is O(1) on the left half. Further, 
C;, contains only the zeros A,,...,A,,. We apply the argument of § 3, taking 
C;, as the path of integration. Let C’, meet the imaginary axis in —ir',, 
ir’. Then r,, r), differ from r, = |A,,| by not more than 3. Instead of L, 


1 iri 
Li, = [flv)dv [ ede. 


0 -il'y 


we now have 


The deviation of C), from circularity will not affect the estimates of § 3. 
It remains to show that —ir), ir}, can be replaced by —ir,, ir, respec- 
tively, with an error which is uniformly o(1). Consider, for example, 
the change from ir}, to ir,. The corresponding change in L‘ is 


1 in 1 : 
[ flv) ae [ ext dz = 25 [ f(x) ME InlE—”) ingen de, 
. P t—t 

where o, = r,—r,), = h(r, +7). 
We may suppose that 6 < 1. Then |o,| < 3}. Since |t—v| < 1, we can 


replace sino, (t—v)/(t—v) by a polynomial 


S (=1)"02"(t—v)" 
ple (2m+1)! 
0 
with an error term less than « if NV > N(e). Hence, by the Riemann- 
This proves that in 





Lebesgue lemma, / is uniformly 0(1) as n> 0. 
Theorem 1 we can put v, = p, and replace r, by |A,,'. 
argument of § 4, we can do the same in Theorem 2. 


In view of the 
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